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ABSTRACT 



Many organizations exist within the Department of Defense 
which exhibit the properties of multi-channel, time-dependent 
queueing systems. One such system is the Army's Acute Minor 
Illness Clinic (AMIC) . Models of this system can be developed 
to determine optimum staffing levels and to upgrade the quali- 
ty of service provided. 

This study developed a time-dependent simulation program 
which was applied to a two queue, multi-channel queueing sys- 
tem. In particular, this model was designed to test and 
evaluate the results of an analytical model of the AMIC. The 
results of both of these models are compared to a second 
simulation program which more closely models the AMIC in order 
to measure the significant differences between it and the two 
other models. 
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I. 



BACKGROUND 



A. AMIC PROGRAM DESCRIPTION 

An analytical model has been developed utilizing an in- 
homogeneous Markov chain and transition matrices to approxi- 
mate solutions of a multi-channel, time-dependent queueing 
system. The system that was modeled is the Acute Minor Ill- 
ness Clinic (AMIC) which is currently in use by the U. S. Army 
at about twenty of their hospitals. The system integrates 
doctors and specially trained physician assistants called 
Amosists into a dual service process. 

Patients are screened as they arrive by a group of Amo- 
sists at a reception desk called the Triage. The nature of 
the complaint is determined and the patients are sent to a 
doctor or Amosist for treatment. A patient may request to be 
seen by a doctor in any event. Patients are treated on a 
first come, first served basis as an Amosist or doctor becomes 
available . 

Amosists treat patients according to set procedures de- 
tailed in technical manuals and occasionally on computer pro- 
grams. These procedures, set up by the staff of physicians, 
indicate the actions an Amosist must take in treating a 
patient, and require the Amosist to seek consultations with a 
physician when the problem is beyond his abilities. This pro- 
cedure helps to eliminate errors in diagnosis and treatment. 
When a consultation is indicated, the Amosist locates the 
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first available physician and together they treat the physi- 
cian. Figure 1 depicts the patient flow of this system. 

B. ANALYTICAL MODEL DEVELOPMENT 

Under the direction of Professor Richard W. Butterworth, 
an analytical model of the stochastic process was developed 
with certain simplifying assumptions of the AMIC. Statistics 
were gathered by David L. VanAsdlen and Leonard 0. Wahlig, 
and an initial model was implemented. The assumptions of 
this model are listed below. 

First, the arrival stream is an inhomogeneous poisson 
process. Log sheets for 75 days were used to determine 
arrival rate function, X(t), for each day of the week. Upon 
examination, it was determined that the arrival function is 
the same for each day of the week with the exception that the 
expected number of arrivals varies. The function is divided 

i 

into 32 time periods with a mean arrival rate for each period. 
Figure 2 is a graph of the input arrival function for the 
model. 

Second, the distribution of services is exponential. The 
Amosist system has a distribution of service times for pati- 
ents that do not see a doctor and a distribution of times for 
patients that are referred. The difference in mean service 
times is a result of the amount of time that an Amosist waits 
with the patient for a doctor to become available and to 
treat the patient. The doctors also have two distributions 
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DIAGRAM OF PATIENT FLOW IN THE REAL CLIN IC 




indicates path of referral patients which preempt 

patients in the physician queue. 



Fif?;ure 1. 
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for service, one for referral patients and one for non-refer- 
rals. The mean of the distribution for physician referral 
patients is considerably shorter than that for patients sent 
directly from the Triage. This is due to the nature of the 
consultation. Mean service times are based on 75 observations 
over a three week period. 

Thirdly, in the model the Amosist and the patient do not 
consult with the doctor. Instead, the queues are kept 
separate and independent by adding the percentage of referral 
patients directly to the physician queue and the Amosist 
queue and computing service times as if the consultation took 
place . 

Let ; 

a = percentage of patients sent directly to the 
physician queue. 

3 = percentage of Amosist patients referred to a 
doctor. 

Y = percentage of patients that are seen by a doctor. 

Then: 



Let: 



Y = a + (1 - a) 3 



E(l) 
E(2) 
E(3) 
E(4) 
Then : 
(1 - a) 3 



mean Amosist non-referral service time 
mean Amosist referral service time 
mean doctor non-referral service time 
mean doctor referral service time 

percentage of referral patients 
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Mean Amosist encounter time = (1-3) E(l) + (3) E(2) 



Mean doctor encounter time; 



(g) E(3) (1 - g)3 (E(4)) 

g+(l-g)3 g+ (l-g)3 



In each case the mean encounter time is the weighted sum of 
the average service times. 

Finally, patients may arrive at the clinic prior to the 
time that the clinic opens. This was modeled under the assump- 
tion that the number of arrivals is distributed poisson, for 
each of the two queues. 

Other parameters for the model include the schedules of 
the Amosists and the physicians, and the volume [italics added] 
or total expected number of patients for the particular day 
under study. For detailed studies of the analytical model, 
see References [1] and [5] . 
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II. SIMULATION MODEL CASE 1 



A. ASSUMPTIONS 

This study developed a simulation program to validate and 
test the significance of the analytical model previously dis- 
cussed. A comparison of results gives an indication of the 
correctness of the model's logic and precision by using the 
same assumptions. 

The simulation uses the same parameters as the analytical 
model. It assumes an inhomogeneous Poisson arrival process 
with three types of patients. These types are used to dis- 
tinguish the number of patients arriving at each location. 
Other than for counting and filing in appropriate queues, 
each patient is considered the same by a server. 

Secondly, patients are served using a mean service rate 
and an exponential distribution for each queue. Patients are 
filed in queues upon their arrival and there is no interaction 
between the queues or the servers. Additional variables are 
added to the physician queue for each referral-type patient 
placed in the Amosist queue. 

Thirdly, for purposes of the model, patients who arrive 
prior to the opening of the clinic begin experiencing waiting 
time when the clinic opens. Figure 3 is a flow diagram which 
depicts patient flow for both the analytical model and the 
simulation model. 
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DIAGRAM OF PATIENT FLOW IN THE 



CASE 1 SIMULATION AND ANALYTICAL MODELS 




= percentage of patients directed to a physician. 
= percentage of patients referred by an Amosist, 






= percentage 
= percentage 
= percentage 



of total patients referred, 
of patients seen by an Amosist only, 
of patients seen by a physician. 



Figure 3* 
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B. MODEL IMPLEMENTATION 



The model is written in Simscript II. 5, a proprietary 
language supported by Consolidated Analysis Centers, Inc., 

12011 San Vicente Blvd. , Los Angeles, California, release 
8D, for implementation on the IBM 360 batch processing system. 
It is an event increment program with first in, first out 
routines. The model generates a daily opening and closing 
of the clinic, beginning and ending empty. As such, each day 
is independent of the others. The simulation is based on a 
particular day of the week for as many iterations as requested. 

Each day the model begins by generating two poisson ran- 
dom deviates. These are the number of patients which arrive 
prior to the opening of the clinic. They are drawn from 
poisson distributions whose means are input parameters to the 
model. A binomial random deviate is then drawn from a dis- 
tribution with the number of trials equal to the number of ' 
arrivals for the Amosist queue and the probability of success 
equal to B/ the percentage of referrals by an Amosist. This 
is the number of Amosist arrivals which are referred to a 
physician. 

Staffing levels change throughout the day as set by input 
parameters to deal with periods of congestion and to allow 
the staff to eat meals. Throughout the day arrivals are gen- 
erated randomly and added to the appropriate queues. 
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The service times are drawn from exponential distribu- 
tions with one mean for each of the queues as discussed in 
the explanation of the analytical model. 

Mean Amosist service time = (1 -3)E(1) + (E(2))3 



Mean doctor service time; 



(g)E(3) ^ (1 -g) 3 (E(4)) 

g + (1 -g) 3 g + (1 -a) 3 



C. ARRIVAL PROCESS 

The arrival stream is generated using the arrival rate 
function supplied to the model. Figure 2 is a graph of the 
function developed from statistics of the Amic. The arrival 
rates for a particular day are computed using the input para- 
meter, volume [italics added] . The volume [italics added] 
is the expected number of patients for a particular day. 

Since part of the daily volume [italics added] arrives prior 
to the opening of the clinic, the arrival rate function does 
not account for all arrivals. Therefore, the arrival rate is 
computed as follows: 

Let : 

V = average number of arrivals per day. 

ave(l) = mean number of arrivals for Amosist queue prior 
to opening time. 

ave(2) = mean number of arrivals for the physician queue 
prior to opening time. 

A(t) = expected number of arrivals in time [0,t]. 

or A(t) = (x)dx 

0 
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where 



A (t) = the arrival rate function shown in Fig. 2. 

Scaling the arrival rate function so as to simulate a given 
average volume of arrivals, v, while keeping the same relative 
shape of the arrival function, gives 

A' (t) = V - ave(l) - ave(2) x A(t) 

T 

where T is the finite horizon (length of the day) . The 
scaled arrival rate function A' (t) is then actually used. 

The actual times of arrival are determined using the cumu- 
lative arrival function. If a set of (x^) of independent ex- 
ponentially distributed variables are drawn with unit mean. 

Let : 



i 




Then S.'s form a homogeneous poisson process. Using the cumu- 
lative arrival function A(t), a relationship exists between 
this homogeneous process and the inhomogeneous process de- 
sired. 

Letting = A ^(S^), the {r^} is one realization of the 
inhomogeneous poisson process and each R^ is the time of an 
arrival. Figure 4 is a graph of the cumulative arrival func- 
tion which is computed by the simulation. 

D. STATISTICS 

Statistics are accumulated for any specific time of the 
day. Because of the nature of the inhomogeneous process, the 
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average for entire days of system statistics, such as system 
size, do not have a useful meaning. The statistics for each 
parameter are the mean and the standard deviation. 



X = E ” X. / N 
i=l ^ 



S 



\/( E X? / N - X^^ 
^ i=l ^ 



The sum and sum of the squares are collected during the simu- 
lation and the means and standard deviations are computed at 
the end. Virtual delay [italics added] is the actual delay 
experienced by each patient during specified time periods. 

The objective is to determine the delay a patient arriving 
at a specific time would experience. In accomplishing this, 
patients who arrive five minutes prior until five minutes af- 
ter a specific time are tagged and their delays are computed 
and averaged. Another statistic, probability that delay is 
greater than an input value called " delay criterion ," [italics 
added] is computed by determining the number of patients in 
the time period that experience delay in excess of that val- 



ue. 



A list of the output statistics follows with definitions: 

Average System Size = Total number of patients in service 

or the queue at a specific time. 

Average Queue Size = Total number of patients in the queue 

at a specific time. 



Utilization Rate = 1 - number of idle servers/number of 

servers on duty. 

(prob that delay=0) = 1.0 if a server is available; 0.0 if 

a server is not available. 
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Virtual Delay 



= delay experienced by patients at a 
specific time. 

(prob that delay > = 1.0 if delay is greater than delay, 

delay . criterion) criterion; 0.0 otherwise. 

E. INPUT TO THE MODEL 

The first part of the simulation includes a routine to 
call input parameters. An example of input parameters is 
appended to this study. Each group of parameters has an op- 
tion word which defines the parameters to be input. The 
routine defaults all parameters to Base Day Monday if input 
is not provided. Option words and their definitions follow. 



Option Word Definition 



Parameter 



ENCOUNTER 



VOLUME . 
TOTAL 



OPENING. 

TIME 



Mean service 1-dimensional array > = 0 

time Encounter (l)=Amosist nonreferral 

mean service time 

Encounter (2)=Amosist referral 
mean service time 

Encounter (3)=doctor nonreferral 
mean service time 

Encounter (4)=doctor referral 
mean service time 



Default values: 

Encounter (1) 

( 2 ) 

(3) 

(4) 



17.96 

23.83 

13.55 

4.34 



Average num- 
ber of pati- 
ents per day 



Integer > 0 

Default value: 
Volume = 156 



24-hour clock real > 0 

time that Default value: 

clinic opens Opening. time = 0800. 
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CLOSING. 

TIME 


24-hour clock 
time that 
clinic closes 


real > 0 
Default value; 
Closing. time = 2400 


PERCENTAGE 


Percentage of 
patients tri- 
aged directly 
to a doctor 


real > 0 
Default value; 
alpha = .55 


REFERRALS 


Percentage of 
patients re- 
ferred by an 
Amosist 


real ^ 0 
Default value; 
Beta = .31 


ITERATIONS 


Number of days 
of simulation 


integer > 0 
Default value; 
Last. day = 1 


PRELOADING 


Average number 
of patients 
arriving prior 
to opening 


1-dimensional array 
Ave(l) = number of Amosist 
arrivals 

Ave(2) = number of doctor 
arrivals 
Default Values; 

Ave . arrival (1) = 7. 

Ave . arrival (2) = 4. 


AMOS 1ST. SCH 


Schedule for 
Amosist staff 


integer > 0 

No.duty.chgs = number of staff 
changes 

Hour(i) = time shift starts 
No. amosists (i) = no of staff 
for that shift 
Hour(l) must = opening. time 
Default values; 

No.duty.chgs = 5 



Hour(l) = 0800 No. ainosists (1) =5 

(2) = 1100 (2)=3 

(3) = 1300 (3)=5 

(4) = 1530 (4)=7 

(5) = 1630 (5)=2 

PHYSICIANS Schedule for integer > 0 

Physician staff No . rosterchgs = number of staff 

changes 

Clock (i) = time shift starts 
No. doctors (i) = number of staff 
for that shift 

Clock (1) must = opening. time 
Default values; 

Clock ( 1 ) =0800 No. doctors (1) =2 
( 2)=1100 ( 2)=1 
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ARRIVAL. FCN 



DELAY . VALUE 



(3) =1200 (3)=2 

(4) =1530 (4)=3 

(5) =1630 (5)=1 

Arrival function integer > 0 

pairs = no of periods de- 
fined 

2 , 1-dimensional arrays 
Time(i) = end of period 
Customer (i) = rate of 
arrival 

alternate time and customer 



Default values: 
Pairs = 32 



time(l) = 810 


Customer (1 ) =22 . 5 


(2) = 820 


(2)=23.5 


(3) = 840 


(3)=24.5 


(4) = 850 


(4)=23.5 


(5) = 900 


(5)=22.5 


(6) = 910 


(6)=21.5 


(7) = 920 


(7)=20.5 


(8) = 930 


(8)=19.6 


(9) = 945 


(9)=18.6 


(10) =1000 


(10)=17.6 


(11) =1015 


(11)=16.6 


(12) =1030 


(12)=14.7 


(13) =1045 


(13)=13.35 


(14) =1100 


(14)=12.05 


(15) =1115 


(15)=10.95 


(16) =1130 


(16)=9.9 


(17) =1145 


(17)=8. 9 


(18) =1200 


(18)=7.9 


(19) =1215 


(19)=7.0 


(20) =1245 


(20)=6.2 


(21) =1300 


(21)=7.2 


(22) =1330 


(22)=13.6 


(23) =1400 


(23)=12.5 


(24) =1430 


(24)=11.1 


(25) =1500 


(25)=9. 85 


(26) =1530 


(26)=8.8 


(27) =1630 


(27)=7.55 


(28) =1730 


(28)=6.3 


(29) =1930 


(29)=5.8 


(30) =2130 


(30)=4.4 


(31) =2330 


(31)=2.8 


(32) =2400 


(32)=. 01 



Time in minutes real > 0 

to compute prob Default value: 

delay is greater Delay . criterion = 15.0 

than delay 

criterion 
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END. OF. DATA End of input. Nothing should follow this 

card. 

The input routine is designed to be as user oriented as 
possible. All input is in free format with a blank space be- 
ing the deliminater for each input word. The input option 
cards and their associated values may be read in whatever 
order desired as long as the end. of. data card is last. 

F. VERIFICATION OF THE PROGRAM 

In the early stages of formulation extensive output was 
used to follow, exactly, the patient flow through the model. 
In its final form the model gives a daily output of the num- 
ber of patients treated in each system by type, the number of 
arrivals prior to opening time for each queue, and the time 
the clinic closed on that day. Provision has been made, if 
the congestion continues throughout the day, to count the 

i 

number of patients that are waiting for service when the 
clinic closes. In the actual clinic these people would be 
sent to the emergency room for treatment. 

The number of referred patients served in both queues was 
always equal and the sum of the referral patients and the 
non-referral patients from both queues always equalled the 
total number of arrivals for the day. 

The means of the poisson deviates and the binomial devi- 
ate from the early arrivals were determined. For patients 
seen only by a doctor the mean was 7.2 with an input mean of 
7 for a 200 day simulation. For the amosist queue the mean 
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was 3.92 with a standard deviation of 1.95 and an input para- 
meter of 4. The expected number of referral patients was 2.17 
and the actual average was 2.24 with a standard deviation of 
1.47. There is no significant difference in the means. 

The total number of arrivals for each queue was also 
summed by type over the simulation to compare to the expected 
values, to insure the correct percentage of patients was 
directed to each queue. 

Ave. num. of doctor's patients = 84.16; expected 

num = 85.8 

Ave. num. of amosist's patients = 48.94; expected 

num = 48.44 

Ave. num of referral patients = 22.36; expected 

num = 21.76 

A test of hypothesis was made on the arrival distribution 
to determine if the number of arrivals for each period was 
from the same distribution. The statistic is Q and is dis- 
tributed x^(n-l) 

32 - nr^^)^ 



where r^ is the mean number of arrivals in the ith period and 

X. is the observed number of arrivals over all iterations. 

1 

2 

At the .025 level of significance the x is 47.0 with 31 de- 
grees of freedom. The test statistic was 37.64, and the hypo- 
thesis was accepted. 

The mean service times were computed using a language 
generated routine, and compared to the expected values. Since 
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the program uses internal exponential generating routines, it 
was only necessary to check the means. 

Amosist mean service time = 19.583; expected value 19.78 

Doctor mean service time = 11.559; expected value 11.687 

The results of the model are tabulated in a later section 
of this study. Source listings and flow charts of this pro- 
gram are appended to the study. 
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III. CASE 2 SIMULATION MODEL 



A. ASSUMPTIONS AND CHANGES TO THE MODEL 

After developing a model of the system using the assump- 
tions of the analytical model, the next step in the study was 
to produce a model whose assumptions more closely duplicated 
those of the clinic. The change required was in the handling 
of referral patients. In this model as in the actual system, 
all patients going to the Amosist have one distribution of 
service times. In the previous model, the extra time for a 
referral patient was time used waiting to see a physician and 
consulting with him. The Amosist is held out of service in 
this model until the consultation is complete and then the 
Amosist can continue treating other patients. In the physi- 
cian queue, referral patients preempt non-referral patients 
and have a separate distribution of service times from the 
patients who are sent directly from the Triage to a doctor. 

All other assumptions of the model remain the same as the 
Case 1 model. Figure 5 is a flow diagram for patients in 
this model. 

B. VERIFICATION 

The input and arrival process for this model was unchanged 
and, therefore, was not retested. The service times averaged 
18.119 minutes for the Amosist queue and 11.419 minutes for 
the physician queue. The expected values were 17.96 minutes 
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DIAGRAM OF PATIFNT FLOW IN THF 



GA'^F 2 SIMULA.TTON VODFL 




indicates path of referral patients and the Amosists 

treatin,e: referral patients. also indicate blocked 

Amosist service channels during a consol tation. 



Figure 5 * 
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for the Amosist queue and 11.687 minutes for the physician 
queue. A source listing for this model is appended to the 
study. 



30 



i 



} 

I 



IV. SIMULATION RESULTS 



A. COMPARISON OF ANALYTICAL AND CASE 1 MODELS 

In order to test the results of the two simulations against 
the analytical model of the AMIC, the means and standard de- 
viations of each of the output parameters were tabulated. The 
results of the analytical model are found in Table I. Table 
II is the results of the Case 1 simulation model. 

Since the simulations and the anlytical model have 
separate distributions and alternate methods of obtaining 
values which are independent, a test can be made that the 
average values are now significantly different. 



Ho: X = A where X = simulation statistic 

A = analytical deterministic value 



For purposes of the test the output parameter of the analyti- 
cal model is assumed to be the true parameter of the system. 
Let; 

N = number of samples in the distribution (2)) 

X = mean value of simulation parameter 
A = mean value of analytical model 
s = standard deviation of X 
T = test statistic 



Then : 



T 



(X - A) 

t: 

S 



X 
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TABLE I. RESULTS OF THE ANALYTICAL MODEL 
AMOS 1ST SYSTEM 
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TABLE I (CONTINUED) 
PHYSICAL SYSTEM 
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RESULTS OF THE CASE 1 SIMULATION MODEL 
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TABLE II.(CONT.) RESULTS OF THE CASE I SIMULATION MODEL 
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Since there are 200 iterations in the simulation, the T sta- 
tistic approaches its limiting distribution which is normal 
(0,1). Tables III and IV are the summary of this comparison. 
All "Undef" entries indicate a standard deviation of 0.0 
indicating that all values of X were equal at that time. 

The two models were compared for the parameters of base 
day Monday. A critical region for accepting or rejecting the 
hypothesis was not chosen since entries are not independent 
of time. Although each day's statistics are independent, the 
statistics taken at different times of the day are serially 
correlated, and cannot be pooled. Therefore, subjective 
judgment based on the normal statistics computed was used to 
judge the fit of the two models. 

Using the above criteria, the analytical model is a very 
good fit in the Amosist system. However, when congestion 
overloads the queues, the analytical model does not conserve, 
all of the probability used in computing the output para- 
meters. As a result, the analytical model underestimates 
the number of patients in the system. Because the physician 
system in this particular set of parameters was highly con- 
gested, the physician system of the analytical model indi- 
cated fewer patients being serviced. The values of both 
models were very close after the congestion diminished. 
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TABLE III. COMPARISON OF THE ANALYTICAL AND SIIVIULATTON MODELS 



NORMAL STATISTICS FOR THE CASE 1 AMOS 1 ST SYSTEM 



System size Queue si 7 .e Utilization Prob(Del= 0 ) Average Delay Time 



- 


.953 


- 


.708 


- 


1.080 


- 


2.295 


- 


. 1.788 


- 


2.287 


— 


2.202 


- 


1.039 


- 


3.056 


- 


1.653 


- 


1.050 


- 


2.098 


mm 


1.320 


- 


1.162 


- 


1.108 


- 


1 . 6 U 6 


- 


1.151 


- 


1.836 




• 997 


- 


.527 


- 


1.405 




.612 


«- 


.202 


- 


1.018 




. 300 




.427 




. 044 




.839 




.737 




.691 




.490 




.448 




.376 


- 


.184 




.088 


- 


.339 




1.16? 


- 


.324 




1.750 


- 


7.575 




.033 




1.37 




5.112 




6.367 


- 


3.096 




4.913 


- 


6.013 


mm 


3. 066 


- 


5.167 


- 


5.325 


- 


3.740 


— 


3.651 


- 


5.957 


- 


2. 020 


- 


1.639 


- 


4.455 


- 


.259 
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2.019 


- 


1.421 


- 


1.801 


.. 


1.384 




Undef 




9.400 


- 


.645 




.274 


- 


.747 




.425 


- 


.429 




.512 


- 


.816 


- 


3.230 


- 


.649 




1.251 




Undef 


- 


1.100 




.583 




Undef 


- 


6.273 


- 


24.392 




Undef 


- 


1.130 




.920 


- 


2.340 




2.746 




2.205 


- 


.973 




3.701 




3.168 




.993 




3.710 




2.910 


- 


.381 




3.185 




3.M7 




Undef 




3.479 





1.107 




.110 


810 




.988 


- 


1.970 


820 




2.045 


- 


.910 


840 




1.318 


- 


3. 166 


850 




.828 


- 


.378 


900 




1.288 


- 


1.208 


910 




2.031 




.920 


920 




1.097 




.222 


930 


- 


. 066 




.699 


945 


- 


1.096 




.003 


1000 


- 


.773 




.356 


1015 


- 


.429 


- 


1.582 


1030 


- 


1.609 




1.406 


1045 


- 


2.8O5 




2.049 


1100 




3.786 


- 


2.830 


1115 




2.067 


- 


5.443 


1130 




4.472 


- 


7.619 


1145 




1.873 


- 


5.100 


1200 




.670 


- 


7.166 


1215 




1.033 


- 


.700 


1245 




15.585 


- 


9.096 


1300 




.225 


- 


.075 


1330 


- 


.537 


- 


• 655 


1400 


- 


.152 




.443 


1430 




Undef 




Undef 


1500 




Undef 




Undef 


1530 


- 


2.176 




1.148 


1630 




.137 




.272 


1730 


- 


1.034 




.963 


1930 


- 


1.696 


- 


.186 


2130 


- 


1.159 




.457 


2330 




Undef 




Undef 


2400 
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TABLE IV. COr.'ipARISON OF THR .ANAT,YTICA L AND S IMULATION fTODET,S 

NORIVIAL STATISTICS FOR THE CASE 1 PHYSICIAN SYSTEM 



System size 


Queue size 


Utilization 


Prob(Del=0) 


Averace Delay 


Time 


13.708 


12.006 




5.069 




10.201 


7.295 


810 


12.263 


10. 527 




11.360 


- 


5. 122 


5.930 


820 


10. 574 


9.541 




1.273 


- 


2.313 


6.273 


840 


10.673 


9.672 




1.455 


- 


5.756 


4. 549 


850 


10. 098 


11.041 




.944 


- 


4. 681 


4.751 


900 


9.872 


10.497 




2.295 


- 


1.182 


5.259 


910 


9.669 


8.923 




.243 


- 


3.665 


5.861 


920 


10.225 


9.474 




.487 




3.486 


4. 870 


930 


9.744 


8.999 




1.035 


- 


.858 


5.688 


945 


9.877 


9.079 




Undef 




3.605 


6.116 


1000 


9.852 


9.156 




1.200 


- 


2.128 


7.019 


1015 


9.681 


8.955 




3.309 


- 


' . 444 


9.278 


1030 


9.432 


8.699 




2.251 


- 


1.748 


7.861 


1045 


7.198 


8.113 




.886 


- 


2.818 


9.489 


1100 


7.486 


7.202 


- 


.424 


- 


. 0^6 


-15.471 


1115 


7.495 


7.189 




.529 




.197 


-16.362 


1130 


7.522 


7.597 




.667 




1.088 


“I9.723 


1145 


7.491 


5.794 




1.325 




.979 


-19. 315 


1200 


7.155 


6.672 




.220 




. Oil 


3.082 


1215 


6.624 


6.013 




3.512 


- 


2.438 


- .137 


1245 


6.078 


5.530 




2.474 


- 


1.515 


2.926 


1300 


6.030 


5.449 




2.124 


mm 


1.962 


1.078 


1330 


6.343 


5.715 




3.503 


- 


1.162 


. 019 


1400 


5. 565 


4.972 




2.730 


- 


2.085 


.757 


1430 


5.598 


4.903 




4.232 




2.226 


- .149 


1500 


4.522 


2.835 




.523 




1.124 


- .487 


1530 


.096 


2.227 


mm 


.257 


- 


5.969 


3.942 


1630 


1 . 062 


2.326 


- 


1.844 




1.720 


- 2.343 


1730 


2.245 


.788 




.271 




1.578 


- .417 


1930 


.727 


.557 


- 


2.245 




1.058 


- 2.079 


2130 


1 . 349 


.909 




.139 




.887 


- 5.231 


2330 


.905 


.874 


- 


.575 




.615 


Undef 


2400 
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B. COMPARISON OF THE CASE 2 SIMULATION AND ANALYTICAL MODELS 



The same procedure outlined for Case 1 was used in com- 
paring the results of the Case 2 model. Table V is the re- 
sults of the Case 2 model and Tables VI and VII are the sum- 
mary of the comparison tests for the models. 

In this case other factors caused differences between the 
two models. In the Case 2 model, Amosists spend more time 
than expected waiting for consultations. Therefore, the two 
models have different distributions for referral service. 

The result was that the Amosist queue was longer than the 
original model and the analytical model results were not as 
close as previously seen. In the other system, the physician 
queue was smaller than in the Case 1 model. This was the re- 
sult of the phenomenon of shortest processing time. Since all 
of the patients with the shorter processing distribution were 
served first, the size of the queues are expected to be short- 
er. The average time in the system for a patient does not 
change, but because of the preemption fewer people are wait- 
ing for service at any one time. 

C. COMPARISON OF CASE 1 AND CASE 2 MODELS 

As discussed in the previous section, the average service 
time for an Amosist referral patient was longer than expected 
in the Case 1 model or the analytical model. The physician 
system had shorter queues and fewer patients in the system 
at any one time due to the fact that the distribution of 
service times for referral patients was much shorter than the 
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TABLE V. RESULTS OF THE CASE 2 SIMULATION MODEL 
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TABLE V.(CONT.) RESULTS OF THE CASE 2 SIMULATION MODEL 
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TABLF VT. 


COMPARISON 


OF THF ANALYTICAL AND SIMULATION MODELS 




NORMAL STATISTICS FOR THF r 


:ASE 


2 AMOSIST 


SYSTEM 




System size 


Queue size 


Utilization 


Prob(Del=0) Averaa:e Delay 


Time 


2.362 


1.979 


2.498 




1.538 


2.856 


810 


3. 048 


2.536 


3.171 


- 


2.687 


4.152 


820 


3.561 
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3.619 


- 


3.391 
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840 


2.712 


2.375 


2.508 


mm 
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3.480 
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2.767 


- 


1.876 


3.622 


900 


3.230 


2.627 
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- 


2.706 


2.153 


910 


2.470 


1.631 


3.276 


- 


1.590 


- .126 


920 


2.631 


1 . 656 


3.359 


•E3 


1.931 


1.078 


930 


2.207 


1.392 


2.676 


- 


2.639 


1.617 


945 


2.411 


1.204 


3.051 


- 


3.388 


3.682 


1000 


3.023 


2.368 


2.841 


- 


3.513 


4.752 


1015 


3.651 


2.453 


3.887 


- 


3.333 


2.964 


1030 


3.^79 


2.447 


3.486 


- 


2.561 


2.644 


1045 


- 2.864 


2.399 


2.408 


- 


4.461 


3.312 


1100 


- . 514 


.910 


— 2 « 586 




2.477 


.644 


1115 


.658 


.670 


.408 




.534 


.048 


1130 


1.532 


1.220 


1.369 


- 


1.876 


.938 


1145 


2,677 


2.275 


2.293 


- 


2.187 


1.356 


1200 


1.117 


1.030 


.843 


- 


.277 


- 1.292 


1215 


1.545 


.672 


1.670 


- 


1.149 


- .480 


1245 


- .237 


- 2.391 


- 8.769 




10.011 


-24.455 


1300 


2.649 


.916 


2.687 


- 


2.612 


1.535 


1330 


.530 


- .429 


. 626 




.410 


- 1.220 


1400 


2.888 


.629 


2.950 


- 


1.319 


.129 


1430 


1.295 


1.179 


1.227 


- 


2.077 


.758 


1500 


1.795 


Undef 


- 4.160 




Undef 


- 1.160 


1530 


-22.622 


Under 


- .284 


- 


2. 398 


0. 0 


1630 


1.995 


.015 


2.608 


- 


1.439 


.763 


1730 


3.895 


1.336 


5.241 


- 


2.822 


1 . 069 


1930 


4.533 


1.970 


5.225 


- 


2.958 


1.173 


2130 


4.828 


2.132 


5. 049 


- 


2.399 


.236 


2330 


4.239 


.717 


4.446 




.120 


Undef 


2400 
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tablp: vit. 



C OMPA R ISON OF THE ANALYTICAL AND STMULATTON MODELS 



NORMAL STATISTICS FOR THE CASK 2 PHYSICIAN SYSTEM 



System size Queue size Utilization 



6 . 6?4 


4.823 


2.048 


6.112 


4 . 506 


2.048 


5.479 


4.271 


- .781 


5 . 1^5 


3.965 


1.455 


5.740 


4, 664 


1.080 


5.609 


4.523 


6 . 021 


5.699 


4. 645 


5.536 


5.163 


4.230 


1.901 


5.031 


4 . 088 


1.486 


5.209 


4.321 


1.135 


5. 656 


4.801 


1.862 


5.405 


4.507 


7.839 


4.925 


4.158 


2.296 


3.127 


4.258 


1.052 


3.034 


2.678 


- .424 


3.309 


2.939 


- .151 


3.299 


3.331 


1.673 


3.212 


1.075 


- 1.405 


2.693 


2.292 


- .951 


2.930 


2.270 


1.024 


2.606 


1.904 


1 . 661 


3.093 


2.364 


1.252 


3.292 


2.479 


2.555 


3.404 


2.623 


2.946 


3.278 


2.693 


.714 


3.325 


1.446 


- .972 


- 1.100 


1.351 


1.454 


.114 


- .206 


- 2.284 


1.786 


1.287 


- .427 


.640 


. 088 


.590 


.862 


.150 


.139 


.464 


- .172 


. 536 



Prob(Del= 0 ) 


Average Delay 


Time 




.922 


7.422 


810 




1 . 644 


6.953 


820 


- 


.607 


5.834 


840 


- 


2.191 


6.620 


850 


- 


1.404 


6.132 


900 


- 


4 . 084 


5.486 


910 


- 


3.665 


7.186 


920 


- 


.869 


6.526 


930 


- 


1.771 


5.579 


945 


- 


.939 


7.100 


1000 


- 


2.128 


8.575 


1015 


- 


4 . 980 


14.987 


1030 


- 


1.748 


9.801 


1045 


- 


2.076 


8.575 


1100 




. 465 


- 14 . 066 


1115 




.197 


-19.519 


1130 


- 


1.593 


-20.656 


1145 




1.862 


-19.765 


1200 




. 817 


2.826 


1215 


- 


.802 


.796 


1245 


- 


. 231 


.984 


1300 


- 


.669 


1.018 


1380 


- 


2.230 


1.633 


1400 


- 


2.085 


.327 


1430 


- 


.116 


- 1.282 


1500 




1.703 


.490 


1530 


- 


7.160 


4.372 


1630 




2.314 


- 2.672 


1730 




.458 


- .123 


1930 




.620 


- 1.768 


2130 


- 


.111 


- 3.957 


2330 




1.071 


Undef 


2400 
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distribution for other patients. The result of servicing 
these shorter processing times first, causes the shorter 
queues. However, the output statistics of the two models were 
quite close, indicating that the Case 1 model is a good model 
of the real system. The fact that the output deterministic 
values for the analytical model were so close to the other 
two models is an indication that this analytical model is a 
good approach to solving the problem. 
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V. CASE 3 SIMULATION MODEL 



A. ASSUMPTIONS 

The weakest of the assumptions in all of the models was 
the assumption of exponential service times distributions. A 
third simulation was conducted to determine the sensitivity 
of the model to the service distributions. A degenerate dis- 
tribution was used for these service times. The constant 
value was equal to the means of the exponential service dis- 
tributions used in the other models. All other assumptions 
remained the same as in the Case 1 model. Table VIII is a 
summary of the results of the simulation. 

B. RESULTS 

Based on this one variation to the base case, the exponen 
tial assumption is apparently robust. Due to the randomness 
inherent in simulations the queues were larger than the previ 
ous models in the earlier part of the day and smaller in the 
later part of the day. There is not sufficient difference to 
indicate that the exponential distribution assumption is in- 
valid. The means of the service times is critical, however, 
as was indicated in the previous section. 
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TABLE VIII. RESULTS OF THE CASE 3 SIMULATION MODEL 
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TABLE VIII.(CONT.) RESULTS OF THE CASE 3 SIMULATION MODEL 
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VI. CONCLUSIONS 



A. COMPARISON OF THE ANALYTICAL MODEL TO THE SIMULATION MODELS 

The purpose of this study was to devise a simulation for 
multi-purpose queueing systems which could be applied to many 
varied systems. Second, the purpose was to determine the pre- 
cision and correctness of an analytical model which had been 
applied to such a system. A first case study, duplicating the 
assumptions of the model indicated that, although some weak- 
nesses existed in the analytical model, the application was a 
good approximation of the AMIC Clinic. The model has a finite 
queue and, therefore, underestimates the system at very high 
levels of congestion. Also, in the real system, the preemp- 
tion of patients by patients whose average service times is 
so much less causes smaller queues than a first in, first out 
system. 

The second simulation program revealed the difference be- 
tween the real system and the analytical model in that the 
real world referral patients preempt patients waiting to see 
physicians, and the Amosists waiting with the patient cannot 
service other patients. The distribution of service times 
for these Amosist referral patients has a larger mean than 
originally believed, causing Amosist statistics for the 
second program to be larger. This program also indicated 
smaller physician statistics due to the preemption. 
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B. SENSITIVITY OF THE MODEL 



A third simulation model indicated that a service distri- 
bution assumption of exponentially distributed service times 
was apparently robust. However, the system is very sensitive 
to the means of the service distributions. In all the analy- 
tical model appears to be a good fit, and useful in providing 
the user with information to upgrade service and optimize 
staffing levels. 
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VI. GENERAL NATURE OF THE SIMULATION MODEL 



The simulation programs were developed in as general a 
way as possible to make them applicable to many varied sys- 
tems. By manipulating the input parameters, any queueing sys 
tern with up to two queues and some limitations on interaction 
is capable of being modeled by these programs. For example, 
a system with only one queue is possible by setting the per- 
centage input parameter to 1.0, thus sending all arrivals to 
one queue. Several assumptions remain. Services are distri- 
buted exponentially, and the arrivals make up a poisson pro- 
cess. The preloading or allowing of arrivals prior to open- 
ing may be surpressed by setting the average arrival prior to 
opening to 0.0. The arrival process may be homogeneous or in 
homogeneous depending on the functional values which are in- 
put. Supplying one arrival rate for the entire day produces 
a homogeneous poisson process, v/ith a mean arrival rate equal 
to the input parameter. 
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CASE 1 Slfy^ULATION FLOW CHART 



MAIN PROGRAM 
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INPUT ROUTINE 



/initial 
/ conditions 



Print 

input 

para- 

meters 





5 ?. 



1 . 




SIMULATION TIMING ROUTINE 
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EVENT OPEN THE CLINIC 



Schedule 
event close 
the clinic 



Determine 

early 

arrivals 



file 

patients ir. 

proper 

queues 



Schedule 

event 

arrival 



Schedule 
first 
Amosist 
shift event 



Schedule 
first 
physician 
shift event 



Schedule 

event 

utilizatior 



Ret. 
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CHANGE AWOSTST SCHEDULE EVENT 



CHANGE PHYSICIAN SCHEDULE EVENT 




Create 

additional 

servers 



Create 
servers for 




File server 
in server ' 
queue 


first shift 


p 


Destroy 




x" — ^ 


excess 


/ 




servers 


( 


vP 






0 — 


File server 
in server 
queue 






Schedule 
a sei'vice, 




t>-^ 


N^n queueP^“^ ^ 


remove 
server from 
^Queucte 




Schedule this 
event for 
next change 






Ret . 



The flov/ charts for thses two events are exactlv the same. 
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EVENT ARRIVAL AND ARRIVAL FUNCTION 




Arrival Function 




Drav/ 

random 

number 

> 

Compute 

poisson 

arrival 

‘ 1 

transform 

to 

inhomos’eneouj. 

arrival 

return to 
arrival 
event with 
arrival time 
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EVENT AMOSTST APPOINTMENT 



EVENT PHYSICIAN APPOINTMENT 




The Amosist and physician service events are identical. 

\ 
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EVENT UTILIZATION 



/Set 
/ period 
( number 





Determine 




system size 


) — ^ — 


for this 




period 



Prob delav 
equals 0 i 
1.0 



Determine 
physician 
system size 
this perioc. 



I 



Determine 
physician 
queue size 
this period 



Commute io 
of ' 

physicians 
in service 













prob delay 
equals 0 

is 0.0 






Determine 
queue size 
for this 
period 



Comoute 
io of 

Amosists in 
service 




Prob dealy 
equals 0 
is 1.0 



X 



Reschedule 
this event 
for next 
time period 



^^^et. 
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EVENT CLOSE THE CLINIC 



PRINT 7 
# of / 

patients 
.served / 








Empty 

queue: 


all 








"patient^ 
not sepv^ed^ 



Count # oi’ 
patients nc 
served 



/Print max 

patien t/j ^ 
Dhysician I 
queue 



Print # 
of 

' patients 
lot servec 




reset 

initial 

conditions 



rint ave, 
service , 
ave ^ of 
patients 

I Z~~1 

Schedule 
event end 
of 

simulation 



Advance to 
next day 
( increment) 



ScViedule 
event 
open the 
clinic 



> 



\/ 



Ret, 
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ill 



EVENT END OF SIMULATION 



'Print 
ff of 
/iterations 



Calculate 
means, std 
deviations 



'Print 
Amosist 
/system 
statistic/ 



Print 7 
^physician/ 
system / 
/s tatistica 

1 ■ 

Sprint / 
ictual # J 
'of arrivals 
)er period 

— r 

calculate 
actual 
arrival 
rate 



Plot 
actual 
arrival 
rate vs 
■ ti - m^ - 
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PREAMBLE 

LAST COLUMN I S 60 
TEMPORARY ENTITIES 

EVERY AMOSIST MAY BELONG TO AN AV A I L . NEXT . AMOS AND HAS 
AN OFFICE 

EVERY PHYSICIAN MAY BELONG TO A NEXT . AVAI L. DOCTOR AND HAS 
A CUBICAL 

EVERY PATIENT MAY BELONG TO AN AMOS. QUEUE, A PHYS. QUEUE, 
HAS A TYPE AND A T I ME . OF . A RR I V AL 



EVENT 



EVERY 

EVERY 

EVERY 

EVERY 



NOTICES INCLUDE END . OF . S I MU L AT I ON , CLOSE 
OPEN. CLINIC, TRl AGE AND UTILIZATION 
AMOS. APPOINT HAS A NUMBER 
PHYS. APPOINT HAS A LISTING 
CHG. IN. MEDICS HAS A SUBSCRIPT 
CHG.NUM.PHYS HAS AN INDEX 



CLINIC, 



THE SY 
DEFI 

PRI ORI 

DEFI 



DEFI 

DEFI 



DEFINE 



DEFI 

DEFI 

DEFI 



STEM Owns A PHYS. QUEUE, AN AMOS. QUEUE, A 
NEXT .AVAI L. DOCTOR AND AN AVA I L .NEXT . AMOS 
AVAI L. NEXT. A'-IOS, NEXT . AVAI L. DOCTOR, PHYS 
AMOS. QUEUE AS FIFO SETS WITHOUT FF,FB,FA 
ROUT INES 

ORDER IS END.OF.SI MULATION,CHG.IN.MEOICS 
CHG.NUM.PHYS, AMOS. APPOINT, PHYS. APPOINT, U 
CLOSE. CLINIC AND TRIAGE 

NUMBER, LIST I NG, OFFICE, CUBICAL , PAIRS, INDE 
SUBSCRIPT, NO NREF ,CONSOL, NONCONSOL, DAY , NO 
TYPE, VOLUME, RATE. NO, NO. DUTY. CHGS AND NO. 
AS I NTEGER VARIABLES 

NUM, NO. DOCTORS, NO. AMOS I STS AS INTEGER,!- 
ALPHA, BETA, CLOSING .TIME , TIME .OF.ARRI VAL, 
HE AN. AMOS. SERV I CE,MEAN.PHYS. SERVICE, INTE 
DELAY. CRITERION AS REAL VARIABLES 
S, U, Ti ME. TC, DOCTOR.T IME. TC , FUNCTION, CUST 
AVE. ARRIVAL, PLOT AND R A TE . F UNCT I ON AS 1- 
PERCENT AS A REAL FUNCTION 
ARRIVAL AS A REAL FUNCTION 
DPLTP AS A FORTRAN ROUTINE 



NE 



TY 



NE 



NE 

NE 



NE 

NE 

NE 



.QUEUE AND 
AND RL 



TILIZATION, 

X, REFERAL, 

, LAST. DAY, 
ROSTER. CHGS 

DIM ARRAYS 
RECORD, 

GRAL AND 

OMER, TI ME, 
DIM ARRAYS 



TALLY AVE. SERVICE. A AS THE MEAN OF A. SERVICE 
TALLY AV. SERVICE. P AS THE MEAN OF P. SERVICE 

DEFINE A. SERVICE AND P. SERVICE AS REAL VARIABLES 
TALLY ACASE AS THE MEAN, SCASE AS THE STD. DEV AND MCASE 
THE MAXIMUM OF CASES 
DEFINE CASES AS AM INTEGER VARIABLE 
TALLY MAXI AS THE MAXIMUM OF N. PHYS. QUEUE 
TALLY MEANl AS THE MEAN AND STD AS THE STD. DEV OF MAXD 
DEFINE MAXD AS A REAL VARIABLE 
DEFINE A. NO. DE L AY, P .NO. DELAY , A. DEL AY, P. DEL AY , 

A. ME AN .DELAY, P .MEAN.DELAY ,P .UTIL, A.SYS. S 
P. SYS. SIZE, A. AV .QUEUE AND P.AV. QUEUE AS 
ARRAYS 

DEFINE MEAN AND STAT AS 2-OIM ARRAYS 
DEFINE COUNTl, C0UNT2 AS I NTEG ER , 1 -D IM ARRAYS 

BEFORE REMOVING FROM AMOS. QUEUE, CALL ACALC 
BEFORE REMOVING FROM PHYS. QUEUE, CALL PCALC 
END 



A. UTIL, 

IZE, 

1-OIM 



CASE 1 SIMULATION MODEL 
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i 



I 



I 



MAI N 

CALL INPUT 
RELEASE INPUT 
ADD I TO DAY 
SKIP 1 OUTPUT LINE 
PRINT 1 LINE THUS 

SIMULATION BEGINS 
SCHEDULE AN OPEN. CLINIC NOW 
START NEW PAGE 
START SIMULATION 
END 



EVENT OPEN. CLINIC SAVING THE EVENT NOTICE 
DEFINE I AS AN INTEGER VARIABLE 
DEFINE X AS A 1-DIM ARRAY 
RESERVE X(«) AS 3 

SCHEDULE A CLOSE. CLINIC AT CLOSING. TIME 
IF AVE. ARRIVAL! 1 )=0 JUMP AHEAD ELSE 

LET NUM(l) =POISSON. FIAVE. ARRIVAL! 1) ♦!) 

IF NUM(l) <= 0 JUMP AHEAD ELSE 

LET NUM(3)=3INOHIAL.F(NUM( 1) , BETA, li 
LET NUMIl )= (-NUM( 3) ) «-NUM( 1 ) 

HERE 

IF AVE.ARRIVAL(2)=0 JUMP AHEAD ELSE 

LET NJM(2) = POISSON. FIAVE. ARRIVALI2) ,1) 

HERE 

FOR 1=1 TO 3 DO 

LET XI I )=REAL.F( NUM( I ) ) 

LOOP 

FOR 1 = 1 TO NUM( I ) <-NUM( 2 ) <-NUM(3) , DO 

LET L0A0 = 1.0/{X(1 )fX( 2)4-X(3)) 

CREATE A PATIENT 

LET T=PERCENT 

IF T<=X( UNLOAD 

LET TYPE ( PATIENT )=i 

LET X( i)=xm-i 

FILE PATIENT IN AMOS. QUEUE 

CYCLE 

ELSE 

IF T<=(X(1) <-X( 2) )*LOAD 

LET TYPE! PATIENT)=2 

LET X(2)=X(2)-1 

FILE PATIENT IN PHYS. QUEUE 

CYCLE 

ELSE 

LET TYPE(PATIENT)=3 

LET XI3) = X(3)-l 

FILE PATIENT IN AMOS. QUEUE 

FILE PATIENT IN PHYS. QUEUE 

LOOP 



SCHEDULE A TRIAGE AT ARRIVAL 

SCHEDULE A UTILIZATION AT FUNCTION! 1) 

SCHEDULE A CH G. I N . MED I C S NOW 

SCHEDULE A CHG.NUM.PHYS NOW 

LET INDEXICHG.NUM. PHYS ) = 1 

LET SUBSCRIPT !CHG. IN. MEDICS ) = 1 

DESTROY OPEN. CLINIC 

RETURN 

END 
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EVENT CHG. IN. MEOICS SAVING THE EVENT NOTICE 
NORMALLY MODE IS INTEGER 
LET I = SUBSCRIPT (CHG. IN. MEOICS I 

IF TIME.V=0.00 

FOR J = 1 TO NO. AMOSISTSI 1) , DO 

CREATE AN AMOSIST 

LET CUBICAL (AMOSIST)=J 

FILE AMOSIST IN AV A I L. NEXT . AMOS 

LOOP 

GO TO APPOINTMENT 

ELSE 

LET MEDICS=NO. AMOSISTSI I ) -NO. AMOSISTSI I-l ) 

IF MEOICS< 0 

FOR EACH AMOSIST OF AV AI L . NEX T . AMOS » DO 

IF CUBICAL (AMOSIST) >NC . AMOS I S T S ( I) 

REMOVE AMOSIST FROM A VA I L . NEXT . AMOS 
DESTROY AMOSIST 

ALWAYS 

LOOP 

JUMP AHEAD 
OTHERWISE 

FOR J = MO. AMOS ISTS { I-l ) + l TO NO . AMOS I STS ( I ), DO 

CREATE AN AMOSIST 

LET CUSICAL(AMOilST) = J 

FILE AMOSIST IN A VA I L . NEXT . AMOS 

LOOP 

' APPCI NTKENT' 

IF THE AMOS. QUEUE IS NOT EMPTY 

LET AMOSIST = F . A V A I L . NEXT . AMOS 

REMOVE AMOSIST FROM AVA I L . flEXT . AMOS 

SCHEDULE AN AMOS. APPOINT NOW 

LET NUMBER! AMOS. APPOINT ) =CU BIC AL( AMOS I ST) 

DESTROY AMOSIST 

ALWAYS 

HERE 

ADD 1 TO I 

LET SUBSCRIPT! CHG. I N. MEOICS ) = I 
IF I <=NO.DUTY.CHGS 

THEN IF TIME. TC( IKCLOSING. TIME 

SCHEDULE THIS CHG . I N. MED I C S AT TIME. TCI I) 

ALWAYS 

RETURN 

END 



63 



EVENT CHG.NUM.PHYS SAVING THE EVENT NOTICE 
NORMALLY MODE IS INTEGER 
LET I = INDEX(CHG.NUM.PHYS) 

IF TIME.V=0.00 

FOR J = I TO NO.DOCTORSt 1 ) ,00 

CREATE A PHYSICIAN 

LET OFFICE (PHYSICIAN) = J 

FILE PHYSICIAN IN NE X T . A VA I L . DOCTOR 

LOOP 

GO TO APPOINTMENT 

ELSE 

LET OOC=NO. DOCTORS! I ) -NO . DOCTORS ( I-l) 

IF DOC < 0 

FOR EACH PHYSICIAN OF NE X T . AV A I L . DOCT OR , DO 
IF OFFIC£( PHYSICI AN) > NO. DOCTORS ( I ) 

REMOVE PHYSICIAN FROM NEX T . A VA I L . DOCT OR 
DESTROY PHYSICIAN 
ALWAYS 
LOOP 

JUMP AHEAD 
OTHERWISE 

FOR J = NO. DOCTORS! I-lH-1 TO NO. DOCTORS ! I) , DO 

CREATE A PHYSICIAN 

LET OFFICE! PHYSICIAN) = J 

FILE PilYSICIAN IN NEX T . A VAI L . DOCT OR 

LOOP 

• APPOI NTMENT* 

IF THE PHYS. QUEUE IS NOT EMPTY 

LET PHYSICIAN = F . NE XT . AV A I L . DOCTOR 
REMOVE PHYSICIAN FROM NEXT . AV AI L . DOCTOR 
SCHEDULE A PHYS. APPOI NT NOW 

LET L I STING! PHYS .APPOINT ) =OFF ICE (PHYSICIAN) 
DESTROY PHYSICIAN 

ALWAYS 

HERE 

ADD 1 TO I 

LET INDEX(CHG.NUM. PHYS)= I 
IF I <= NO. ROSTER. CH3S 

THEN IF DOCTOR. TIME. TC( I XCLCSING. TIME 

SCHEDULE THIS CHG.NUM.PHYS AT DOCTOR .T IME . TC ! I ) 

ALWAYS 

RETURN 

END 
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EVENT TRIAGE SAVING THE EVENT NOTICE 
CREATE A PATIENT 

LET TIME.OF. ARR IVAL( PAT lENT ) = TIME.V 
LET T=PERCENT 
IF T <=ALPHA 

LET TYPE! PATI EMT) =2 
FILE PATIENT IN PHYS. QUEUE 
JUMP AHEAD 
OTHERWISE 

IF T <= ALPHA (1 .0-ALPHA)«( 1.0-BETA) 

LET TYPF( PATIENT) =1 
FILE PATIENT IN AMOS. QUEUE 
JUMP AHEAD 
ELSE- 

LET TYPE! PATI ENT) =3 
FILE PATIENT IN AMOS. QUEUE 
FILE PATIENT IN PHYS. QUEUE 
HERE 

IF TYPE (PATIENT )=2 OR T Y PE ( P AT I EN T ) = 3 
THEN IF NEXT. AVAIL. DOCTOR IS NOT EMPTY 
LET PHYSICIAN = F . N£ XT . A V A I L . DOCTOR 
REMOVE PHYSICIAN FROM NEX T . AV A I L . DOC TOR 
SCHEDULE A PHYS. APPOINT NOW 

LET L 1ST I NG (PHYS. APPOINT ) =OFF ICE ( PHYSICIAN) 
DESTROY PHYSICIAN 
ALWAYS 

IF TYPE(PATIENT)=1 DR T YPE ( P A T I EN T ) =3 

THEN IF THE AV A I L . NE X T . AMOS IS NOT EMPTY 
LET AMOSIST = F. AVA I L . NEX T . AMOS 
REMOVE AMOSIST FROM A V A IL . NEX T . AMOS 
SCHEDULE AN AMOS. APPOINT NOW 
LET NUMBER ( AMOS. APPO I NT ) = CUB I C AL ( AMOS I ST ) 
DESTROY AMOSIST 
ALWAYS 

LET NEXT = ARRIVAL 
IF NEXT < CLOSING. TIME 

SCHEDULE THIS TRIAGE AT NEXT 
RETURN 
OTHERWI SE 
DESTROY TRIAGE 
RETURN 
END 
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EVENT AMOS. APPOINT SAVING THE EVENT NOTICE 
DEFINE I AS AM INTEGER VARIABLE 
LET I=NUMBEPx( AMOS. APPOINT ) 

IF I > NO. AMOSI STS( SU3SCRIPT(CHG. IN.MEDICSI-l 1 
DESTROY AMOS. APPOINT 
JUMP AHEAD 

ELSE 

IF THE AMOS. QUEUE IS EMPTY 

CREATE AN AMOSIST 

LET CUBICAL! AMOS 1ST) =NUMBER( AM OS. APPOINT) 
FILE AMOSIST IN AV A I L .NEXT . AMO S 
DESTROY AMOS. APPOINT 
RETURN 

F. AMOS .QUEUE 

= - ME A N. AMOS. SERVICE* LOG. E.F( PERCENT) 



ELSE 

LET PATIENT = 
LET A. SERVICE 



IF TYPE(PATIENT)=1 

ADD 1 TO NONREF 

ALWAYS 

IF TYPE (PATIENT>=3 

ADD 1 TO REFERAL 

ALWAYS 

REMOVE PATIENT FROM AMOS. QUEUE 



IF M.PHYS.QUEUE(PATIENT) = 0 
ALWAYS 



DESTROY PATIENT 



IF A. SERVICE +TIME.V < CLOSING. TIME 

SCHEDULE THIS AMOS. APPOINT AT A . SERV ICE 4-T I ME . V 

REGARDLESS 

HERE 

IF THE AMOS. QUEUE IS NOT EMPTY 

FOR 1=1 TO N. AMOS. QUEUE 

WHILE AVAIL. NEXT. AMOS IS NOT EMPTY, DO 

LET AMOSIST = F . A VA I L . NEXT . AMO S 

REMOVE AMOSIST FROM AVA I L . NEXT .AMOS 

SCHEDULE AN AMOS. APPOINT NOW 

LET NUMBER ( AMO S. APPO I NT ) =CUBICAL( AMOSIST) 

DESTROY AMOSIST 

LOOP 



ALWAYS 
RETURN END 
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EVENT PHYS. APPOINT SAVING THE EVENT NOTICE 
DEFINE J AS AN INTEGEt^ VARIABLE 
LET J = LISTING(PHYS. APPOINT) 

IF J > NO. DOCTORS! INDEX(CHG.NUM.PHYS)-l) 

DESTROY PHYS. APPOINT 
JUhP AHEAD 

ELSE 

IF THE PHYS. QUEUE IS EMPTY 

CREATE A PHYSICIAN 

LET OFFICE ( P HY S I C I A N ) =L I ST I NIG ( PHYS . APP 0 I NT ) 
FILE PHYSICIAN IN N EXT . AV A I L . DOC TOR 
DESTROY PHYS. APPOINT 
RETURN 

ELSE 

LET PATIENT = F. PHYS. QUEUE 

LET P. SERVICE = -ME AN . P HY S . SERVI CE*LOG. E . F ( PERCENT ) 

IF TYPE(PATIENT)=2 

ADD 1 TO NONCONSOL 

ALWAYS 

IF TYPE(PATIENT)= 3 

ADD 1 TO CONSOL 

ALWAYS 

REMOVE PATIENT FROM PHYS. QUEUE 

IF M. AMOS. QUEUE! PATIENT) = 0 
DESTROY PATIENT 

REGARDLESS 

IF P. SERVICE + TIME.V < CLOSING. TIME 

SCHEDULE THIS PHYS. APPOINT AT P . SERVICE «-T I ME .V 

ALWAYS 

HERE 

IF THE PHYS. QUEUE IS NOT EMPTY 
FOR 1=1 TO N. PHYS. QUEUE 

WHILE NEXT .AVAIL .DOCTOR IS NOT EMPTY, DO 
LET PHYSICIAN = F. NEXT , AVAIL . DOCTOR 
REMOVE PHYSICIAN FRCh NEXT .AVAI L . DOCTOR 
SCHEDULE A PHYS. APPOINT NOW 

LET LIST I NG! PHYS. APPOINT )=OFF ICE {PHYSIC I AN) 

DESTROY PHYSICIAN 

LOOP 

ALWAYS 
RETURN END 



67 



EVENT UTILIZATION SAVING THE EVENT NOTICE 
DEFINE I AS INTEGER, SAVED VARIABLE 
IF TIME.V=FUNCTION( 1 ) LET I =1 REGARDLESS 
LET Y = NO.AMOSI STS( SUBSCRIPTICHG. IN.MEDICSl-1 ) 

LET DUMMY = N. AMOS. QUEUE Y - N . AV A I L . NEXT . A MOS 

ADD DUMMY TO A . SYS . S I Z E ( I ) 

ADO DUMMY DUMMY TO ST ATI 1,1) 

ADD N. AMOS. QUEUE TO A . A V . QUE UE ( I ) 

ADO N. AflOS .QUEUE * N. AMOS. QUEUE TO STAT(3,I) 

LET DUMMY = 1.0- N . AVA I L . NEX T . AMOS/ Y 
ADD DUE MY TO A.UTILI I ) 

ADD DUMMY # DUMMY TO ST AT (5,1) 

IF AVAIL. NEXT. AMOS IS NOT EMPTY 
ADD 1.0 TO A.NO.DELAYd ) 

REGARDLESS 

LET Y = NO. DOCTORS! INDEX (CHG.NUM. PHYS )- 1) 

LET DUMMY = N.PHYS. QUEUE f Y - N . NEXT . AVAI L . DOCTOR 
ADD DUMMY TO P . SYS . S I Z E ( I ) 

ADD DUMMY * DUMMY TO ST AT (2, I) 

ADD N.PHYS. QUEUE TO P . A V . QUEUE ( I ) 

ADD N.PHYS. QUEUE « N.PHYS. QUEUE TO STAT(4,I) 

LET DUMMY = 1.0- M . NEX T . AV A I L . DOCTOR/ Y 
ADD DUMMY TO P.UTILI I » 

ADD DUMMY DUMMY TO STAT(6,I) 

IF NEXT. AVAIL .DOCTOR IS NOT EMPTY 
ADO 1.0 TO P.NO.OELAY( I ) 

REGARDLESS 

LET Y = FUMCTICN(I) 

IF Y <= CLOSING. TIME 
ADD 1 TO I 

SCHEDULE THIS UTILIZATION AT FUNCTION(I) 

RETURN 

ELSE 

DESTROY UTILIZATION 

RETURN 

END 
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EVENT CLOSE. CLINIC SAVING THE EVENT NOTICE 
NORMALLY MODE IS INTEGER 

PRINT 1 LINE WITH NU M ( 1 ) , NUM I 2 ) , NUM ( 3 ) THUS 
PRELOADING FOR QUEUES L TO 3 I S * « * 

PRINT I DOUBLE LINE WITH REFERAL , NCNREF , CONSOL , NONCONSOL 
AND TIME.V THUS 

REFERALS * NONREF « CONSOLS « NONCON * CLOSING 
TIME 

FOR EACH AMOSIST OF AV A I L . NEX T . AMOS , D 0 

REMOVE AMOSIST FROM AV A I L . MEX T . AMOS 

DESTROY AMOSIST 

LOOP 

FOR EACH PHYSICIAN OF NEXT . AVAI L . DOCTOR , DO 

REMOVE PHYSICIAN FROM NEX T . AVA IL . DOCTOR 

DESTROY PHYSICIAN 

LOOP 

IF THE AMOS. QUEUE IS NOT EMPTY 

FOR EACH PATIENT IN AMOS. QUEUE, DO 

REMOVE PATIENT FROM AMOS. QUEUE 
IF M.PHYS.QUeUEIPATIENT) = 0 
DESTROY PATIENT 
ALWAYS 

ADD 1 TO A. PATIENTS. NOT. SEEN 
LOOP 

ALWAYS 

IF THE PHYS. QUEUE IS NOT EMPTY 

FOR EACH PATIENT IN PHYS. QUEUE, DO 

REMOVE PATIENT FROM PHYS. QUEUE 

DESTROY PATIENT 

ADD 1 TO P. PATIENTS. NOT. SEEN 

LOOP 

ALWAYS 

PRINT 1 LINE WITH A . P AT I E NT S . NO T . S E EN AND 
P. PATI ENTS. NOT. SEEN THUS 

LEFT IN AMOSIST CLINIC LEFT IN DOCTORS OFFICE 

PRINT 1 LINE WITH MAXI THUS 

THE MAX NUMBER OF PATIENTS IN THE PHYSICIAN QUEUE IS ^ 

SKIP 1 OUTPUT LINE 
LET MAXD = MAXI 
IF DAY =LAST.OAY 

PRINT 2 LINES WITH AV. SERVICE. P AND AVE. SER V ICE .A THUS 
AVERAGE PHYSICIAN SERVICE AVERAGE AMOSIST SERVICE 

PRINT 1 LINE WITH IMEANl AND STD THUS 
THE AVE MAX FDR PHYS. QUEUE IS STD DEV *«*.*«^c 

SKIP 1 OUTPUT LINE 

PRINT 1 DOUBLE LINE WITH ACASE,SCASE AND MCASE THUS 
AN NO OF PATIENTS STD DEV MAXIMUM 

* PER DAY 

SCHEDULE AN END. OF . S I MULAT ION NOW 
RETURN 
ELSE 

LET TIME.V =0.00 LET DAY=DAY<-1 LET NC=0 LET RATE.N0=1 
LET REFERAL=0 LET C0NS0L=0 LET NONCONSOL=0 LET NONREF=0 
LET REC0RD=0 RESET TOTALS OF N. PHYS. QUEUE 
DESTROY CHG. IN. MEDICS DESTROY CHG.NUM.PHYS 
DESTROY CLOSE. CLINIC SCHEDULE AN OPEN. CLINIC NOW 

RETURN 
END 
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EVENT END. OF. SIMULATION 

DEFINE MODCUR AS AN INTEGER VARIABLE 
PRINT 1 LINE WITH DAY THUS 
END. OF .SIMULATION DAY IS 

FOR J = 1 TO 8 ALSU FOR I = 1 TO PAIRS, DO 

IF J = 7 OR J = 8 LET STATU, I) = MEANU,!) REGARDLESS 
LET MEANU, I) = MEAN ( J , I j / L AS T . D AY 
LET STATU, I) = ST ATU , I ) /LAST .DAY-MEANU, I )’!‘MEANU, I) 
LET STATU, IJ = SORT. F ( STAT U , I H 
LOOP 

FOR I = 1 TO PAIRS, DO 

IF COUNTim = 0.0 JUMP AHEAD ELSE 
LET STAT (9,1) = MEAN(9, I ) 

LET ME AN (9, I) = M E AN ( 9 , I) /C OUNT 1 ( I ) 

LET STAT(9,I)=SQRT.F(STAT(9,I)/C0UNTim- MEAN(9,D* 
MEAN(9, I ) J 

LET MEANdl,!) = ME AN ( I 1 , I ) / COUNT 1 U ) 

LET STAT(11,I)=SQRT.F(STAT(11,I )/ COUNT 1( I ) -MEANdl, I )« 
MEANdl, I d 

HERE 

IF C0UNT2(I) = 0.0 CYCLE ELSE 
LET STAT(10,I ) = MEAN( 10,1 ) 

LET M£AN(10,d = MEANdO,!) /COUNTZd) 

LET STAT( 10, I)=SQRT.F(STAT( 10,1 )/CCUNT2( n-MEANdO, I) * 
MEAN( 10,1)) 

LET MEAN(12,I) = M£AN(12, I ) /COUNT Z( I ) 

LET STATd2,I)=SQRT.F(STAT(12,I)/COUNT2( I)-MEAN(12, I)* 
MEAN (12,1 ) ) 

LOOP 

START NEW PAGE 
PRINT 1 LINE THUS 

AMOSIST SYSTEM 

SKIP 1 OUTPUT LINE 

PRINT 1 double line WITH DEL AY . C RI T£ R I ON THUS 

SYS SIZE STD DEV AVE QUEUE STD UEV UTILIZ STD DEV 
P(C=0) STD DEV P(D> AVE DELAY STD DEV NO TIME 

SKIP 1 OUTPUT LINE 
FOR I = 1 TO PAIRS, DO 

PRINT 1 DOUBLE LINE WITH ME AN ( 1 , I ) , ST AT ( 1 , I) , ME AN ( 3 , I ) , 
STAT(3, I ) ,MEAN( 5, 1 ) , STAT(8, I ) ,MEAN(7, I ) , STAT( 7, d , 
MEAN (9, I ) ,MEAN( 11, I), ST AT (11, I) ,C0UNT1( I ) AND 
TIMEd) THUS 

IF FRAC.F( 1/5 .0) =0.0 SKIP 1 OUTPUT LINE ALWAYS 
LOOP 

START NEW PAGE 
PRINT 1 LINE THUS 

PHYSICIAN SYSTEM 

SKIP 1 OUTPUT LINE 

PRINT 1 DOUBLE LINE WITH DEL AY . C R I TE R I ON THUS 

SYS SIZE STD DEV AVE QUEUE STD DEV UTILIZ STD DEV 
P(C=0) STD DEV P(D> =^.) AVE DELAY STD DEV NO TIME 
SKIP 1 OUTPUT LINE 
FOR I =1 TO PAIRS, DO 

PRINT 1 DOUBLE LINE WITH ME AN( 2 , 1 ) , STAT ( 2 , 1 ) , MEAN ( 4 , I J , 
STAT(4, I ) ,MEAN(t), I ) ,STAT (6, I ) ,MEAN(8 , I ) , STAT( 8, I ) , 
MEAN! ID, I ) ,MEAN(12, I) , STAT(12 ,9) ,C0UNT2( I ) AND 
TIMEd) THUS 

IF FRAC.F( 1/5.0) =0.0 SKIP 1 OUTPUT LINE ALWAYS 
LOOP 

START NEW PAGE 
PRINT 1 LINE THUS 

ACTUAL NUMBER OF ARRIVALS PER PERIOD 
SKIP 1 OUTPUT LINE 

FOR 1= 1 TO PAIRS, WRITE PLOT(I) AS D (15,4) 

FOR 1= 1 TO PAIRS- 1, DO 

IF I = 1 LET PLOT! I )=PLOTd )-60. 0/(FUNCTI oN( 1 J»LAST .DAY) 
ALWAYS 

LET PLOT! H-lJ=PLOT(d-l)«60.0/( (FUNCTION! I + D- 
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FUNCTION( I ) )*LAST.DAY) 

LOOP 

CALL DPLTP( TIME( ,PLOT(*) , PA I R S , MODCUR ) 
STOP END 



ROUTINE PERCENT 
LET X=RANDOM.F(2 ) 
RETURN WITH X 
END 



ROUTINE INPUT 

DEFINE I AS AN INTEGER VARIABLE 

DEFINE ENCOUNTER, HOUR AND CLOCK AS 1-DlM ARRAYS 
DEFINE OPTION AS AN ALPHA, 1-DIM ARRAY 

RESERVE ENCOUNTER(^) ,AVE. ARRIVAL!*) ,NUM(*I AND OPTION!*) 
AS 4 

LET ENCOUNTER!!) = 17.96 LET ENCOUNTER!2) = 23.83 
LET ENCOUNTER!3) = 13.55 LET ENCOUNTER!^) = 4.34 
LET VOLUME = 156 

LET OPENING. TIME = 300. LET CLOSING. TIME = 2400. 

LET ALPHA = .55 LET BETA = .31 LET LAST. DAY = 1 
LET AVE. ARRIVAL! 1 ) = 7. LET AVE . ARR I VAL! 2 ) = 4. 

LET DELAY. CRITERIJN=15.0 
HERE 

IF MODE IS NOT ALPHA 

PRINT 1 LINE WITH READ.V THUS 

EXPECTED OPTION WORD WHILE READING CARD NUMBER * 
STOP 
OTHERWI SE 

IF EFIELD.F-SFI ELD. F > 12 

PRINT 1 LINE WITH READ.V THUS 
WHILE READING CARD NJM3ER * OPTION WORD WAS TOO LONG 
STOP 
OTHERWISE 

READ OPTION! 1 ) tOPTION!2) , OPTION! 3) 

IF OPTION!!) = "ENCO" 

FOR I = I TO 4, DO 

READ ENCOUNTER! I ) 

LOOP 

JUMP BACK 

ELSE 

IF OPTION!I)= "VOLU" 

READ Y 

LET VDLUME= Y 
JUMP BACK 

ELSE 

IF OPTION!1)="OPEN" 

READ OPENING. TIME 
JUMP BACK 

ELSE 

IF OPTION!! )="CLOS” 

READ CLOSING. TIME 
JUMP BACK 

ELSE 

IF OPTI 0N(1 )= "PERC" 

READ ALPHA 
JUMP BACK 

ELSE 

IF OPTION! 1)="REFE" 

READ BETA 
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JUMP BACK 

ELSE 

IF OPTION! 1 ) = "ITER^' 

READ Y 

LET LAST. DAY = Y 
JUMP BACK 

ELSE 

IF OPTION(l)="PREL" 

READ AVE. ARRIVAL! 1 ) , A VE . ARRI VAL ! 2 ) 

JUMP BACK 

ELSE 

IF OPTION! I )="AMOS” 

READ Y 

LET NO.DUTY.CHGS = Y 

RESERVE TIME.TC!*),NO.AMOSISTS!«) AMD CLOCK!*) AS 
NO.DUTY.CHGS 

FOR I = I TO NO.DUTY.CHGS ,DO 

READ CLOCK! I ) ,NO.AMOSISTS! I) 

LOOP 

JUMP BACK 

ELSE 

IF OPTION!! )="PHYS" 

READ Y 

LET NO. ROSTER. CHGS = Y 

RESERVE DOCTOR. TIME. TC!*) , NO. DOCTORS!*) , HOUR!*) AS 
NO. ROSTER .CHGS 

FOR I = 1 TO NO. ROSTER. CHGSt DO 

READ HOUR! I ) ,NO. DOCTORS! I) 

LOOP 

JUMP BACK 

ELSE 

IF OPTICN!l)="ARRI*» 

READ Y 

LET PAIRS = Y 

RESERVE FUNCTION! * ) ,RATE.FUNCT I ONI *) , 
CUSTOMER!* ) jCOUNTl !*) AND COUNT?!*) AS PAIRS 
RESERVE TIME AS PAIRS 

RESERVE STAT!*,*) AND MEAN!*,*) AS 12 BY PAIRS 
FOR I = I TO PAIRS, DO 

READ TIME! I ), CUSTOMER! I ) 

LOOP 

JUMP BACK 

ELSE 

IF OPTI ON!l ) ="DELA" 

READ DELAY.C RITERION 
JUMP BACK 

ELSE 

IF OPTION!!) -,= "END." 

PRINT 1 DOUBLE LINE WITH OPTIONII) THUS 
**** IS NOT THE PROPER BEGINNING TO AN OPTION WORD. SIMULAT 
ION ENDED 

STOP 

OTHERWISE 

PRINT I DOUBLE LINE WITH ENCOUNTER!!), ENCOUNTER!?), 
ENCOUNTER !3) AND ENCOUNTER! -t ) THUS 
ENCOUNTER!!) ***.v« ENCOUNTERI2) ***.** ENC0UNTERI3) ***. 
** ENCCUNTERI4) *.** 

PRINT I DOUBLE LINE WITH VO LUME , 0 PEN I NG . T I ME , CLOS ING.T I ME , 
ALPHA, BETA AND LAST. DAY THUS 

VOLUME * OPENING TIME * CLOSING TIME * PERCENT TO 

DOCTOR .*** PERCENT REFERRED .**" NO OF DAYS * 

LET ME AN. AMOS. SERVICE = ! 1 . 0-6ET A ) *ENCOUNTER ! 1 ) t-B E TA 
*ENC OUNTER! 2) 

LET MEAN.PHYS.SERVICE= ! ALP H A*ENC OUNT E R ! 3) + ! 1 . 0- ALPH A ) * 
ENCOUNTER ! 4) *B ETA) /! ALPHA+! I . 0- ALPHA ) *BETA ) 

PRINT 1 LINE WITH ME AN . AMOS . SER V I C E , M EAN . PHY S . S 5R V I C E THUS 
MEAN SERVICE TIME FOR AMOSIST FOR DOCTOR ***.*** 

SKI? I OUTPUT LINE 
FOR I = 1 TO 2 , DO 

PRINT I LINE WITH I , AVE . A RRI V AL ! I ) THUS 
AVE PRELOAOING OF QUEUE * IS * 

LOOP 
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SKIP 1 OUTPUT LINE 
IF NO.DUTY.CHGS <=0 

LET MQ.OUTY.CHGS = 5 

RESERVE TIME.TCI ) ,NO. AMOSI STS( =«=) AND CLOCK!*) AS 5 
LET NO. AMOSI STS( 1 ) =5 LET NO . AMO S I S TS ( 2 ) =3 LET 

NO. AMOSI STS ( J) =‘3 LET NO . AMO S I S T S ( 4 ) =7 LET 
NO.AMOSi STS(5) =2 LET CLOCK(l) =800. LET CLOCK ( 2 ) = i 1 00 . 
LET CLQCK(3)=1300. LET CLOCK! 4) = 1 530. LET CLOCK ! 5 ) = 1 630. 
REGARDLESS 

IF NO. ROSTER. CHGS <=0 

LET NO. ROSTER. CHGS=5 

RESERVE DOCTOR. TIME. TC!*) , NO. DOCTORS!*) , HOUR!*) AS 5 
LET NO. DOCTORS! i )=2 LET NO . D OC TORS ! 2 ) = 1 LET 
N0.D0CT0RS!:>)=2 LET NO . DOC TORS ! 4 ) =3 LET NO . DOCTOR S ! 5 ) =1 
LET H0UR!1)=80C. LET HOUR ! 2 ) = 1 1 00 . LET HOUR ! 3 ) = I 200 . 

LET H0UR!4) =1530. LET HOUR ! 5 ) = 1 63 0. 

REGARDLESS 

FOR I = 1 TO NO. DUTY. CHGS, 00 

LET Y = CLOCK! I )-OPENING. TIME 
LET Y=Y/iOO.O 

LET TIME.TC!!) = TRUMC . F ! Y ) *60 . 0+FRAC . F ! Y ) *1 00 . 0 
PRINT 1 LINE WITH NO . AMO S I S T S ! I ) , CLOCK ! I ) THUS 
NUMBER OF AMOSI STS ON DUTY *** TIME SHIFT STARTS * 

LOOP 

SKIP 1 OUTPUT LINE 

FOR I = 1 TO NO. ROSTER. CHGS ,00 

LET Y = HOUR! I )-OPENING.TIME 
LET Y= Y/100.0 

LET DOCTOR. TIME. TC! I)=TRUNC.F!Y)*60.0 +FRAC.F!Y)* 
100. 0 

PRINT 1 LINE WITH NO. DOC TORS ! I ) ,HOUR ! I ) THUS 
NUMBER OF DOCTORS ON DUTY * TIME SHIFT STARTS * 

LOOP 

SKIP 1 OUTPUT LINE 
IF PAIRS <=0 

LET PAIRS = 32 

RESERVE FUNCTION! *) ,RATE . FUNCT I ON ! *) ,CUSTOMEP!*) , 
COUNTll*) .C0UNT2!*) AND TIME!*) AS 32 
RESERVE STAT!*,*) AND MEAN!*,*) AS 12 BY 32 
LET TIME!1I=810. LET TIME!2)=820. LET TIME!3)=840. LET 
TIME!4)=850. LET TIME!5)=900. LET TIME!5)=910. LET 
TIME!7)=920. LET TIME!8)=930. LET TIME!9)=945. LET 
TIME! 10) =1000. LET T I ME ! 1 1 ) = 1 01 5 . LET T I ME ! 1 2 ) =1 0 3 0 . 

LET TIME!13)=1045. LET T IME ! 14) = 1 1 00. LET TIME!15)= 
1115. LET TIME! 15)=1130. LET TI ME ! 17 ) = 1 1 45. LET 
TIME! 18) =1200. LET T I ME ! 19 ) = 1 2 1 5 . LET T I ME ! 2 0 ) = 1 24 5 . 

LET TIME !21) =1300. LET T IME ! 22) = I 330. LET TIME!23) = 
1400. LET TIME! 24) =1430. LET T I ME ! 25 ) = 1 50 0 . LET 
TIME!26) =1530. LET T I M E ! 2? ) = 1 630 . LET T I ME ! 2 8 ) = 1 73 0. 

LET TIME ! 29) = 1930 . LET TI ME ! 30 ) =2 130. LET TIME!31) = 
2330. LET TIME!32)=2400. 

LET CUSTOMER! 1) = 22. 5 LET CU STCME R ! 2 ) = 2 3 . 5 LET CUST0MER!3) = 
24.5 LET CUSTOMER! 4)=23. 5 LET CUS TOMER ! 5 ) =22 . 5 LET 
CUST0MER!6)=21.5 LET C US TOM E R ! 7 ) = 20. 5 LET C US T OME R ! 8 ) = 1 9 . 6 
LET CUST0MER!9)=18.6 LET CU STOME R (10 ) = 1 7. 6 LET CUSTOMER(ll) 
=16.6 LET CUSTOMER ! 12) =14. 7 LET C UST OME R ! 1 3 ) = 1 3. 3 5 LET 
CUSTOMER! 14 )= 12. 05 LET C UST OME R ! 1 5 ) = 1 0 . 95 LET CUSTOMER ! 1 6 ) = 
9.9 LET CUST0MER!17)=8.9 LET C USTOME R ! 1 8 ) =7 . 9 LET 
CUSTOMER! 19) =7. LET C U ST OME R ! 20 ) = 6 . 2 LET CUST 0 MER ( 2 1 ) = 7 . 2 
LET CUSTOMER! 22) =13.6 LET CUS70ME R( 23 ) = 12. 5 LET CUST0MER!24) 
=11.1 LET CUST0MER!25)=9.85 LET C USTOMER ( 26 ) = 8 . 8 LET 
CUSTOMER! 27)=7. 55 LET C USTOM ER ! 28 ) =6 . 3 LET CUSTOMER ! 29 ) =5. 8 
LET CUST0MER!30)=4.A LET CU S TOME R ! 3 1 ) = 2 . 8 LET CUS TOMER! 32 ) = 
.01 

REGARDLESS 

FOR 1= 1 TC PAIRS, DO 

LET Y =!TIME! I)-0PENING.TIME)/100.C 

LET FUNCTION! I ) =TRUNC.F!Y)*60.0 *■ FRAC.F!Y) *100.0 
PRINT 1 LINE WITH I , C UST OME R ! 1 ) , T 1 ME ! I ) THUS 
PERIOD NUMBER * AVE ARRIVAL RATE *.**'* AT TIME * 

LOOP 

LET Y = (CLOSING. TIME-OPENING. TIME)/100.0 
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LET CLOSING. TIME = TRJNC.F{Y)*60.0 + F RAC . F ( Y ) 00 . 0 

LET A.SYS.SUEI*) = MEANll,*) LET P . SY S . S I Z E ( ) = MEAN(2,«) 
LET A.AV.QUEUE(>:=) = MEAN(3,*) LET P . A V . QUEUE ( * ) = MEAN(4,*) 
LET A.UTILI^') = MEAN(5,«) LET P.UTIL(*I = MEAN(6,*) 

LET A. NO.D£LAY(« » = MEAN(7t=^> LET P . NO. DEL A Y ( * ) = MEAN(8,*) 
LET A.DELAY(=:») = MEAN(9,=<') LET P.DELAY(=:') = MEAMdO,*) 

LET A.MEAN.DELAY(>!=) = MEAN! 11, «) LET P . MEAN . D EL A Y( « ) = 
MEAN(12,=!‘ J 
IF PAIRS=1 

LET INTEGRAL=CL0SING.TIME*CUST0MER( I ) 

JUMP AHEAD 

ELSE 

LET INTEGRAL = FUNCT I ON ( 1) ^CUSTOM ER ( 1 J 

FOR I =■■ 2 TO PAIRS WHILE FUNCTION! i) < C LOS I NG . T I ME , DO 

LET INTEGRAL = I NT EGR AL + ( FUNCT I ON ( I) -F UMC T I ON ( I - 1 ) ) 
^CUSTOMER! I ) 

LOOP 

LET INTEGRAL = I NTEGR AL + ( CLOS I NG. T I ME-FUNCT I ON (I- 1) ) 
^CUSTOMER! I ) 

HERE 

LET RE MAI MDER=VOLUME-AVE .ARRIVAL! 1 ) -AVE . ARRI V AL ( 2 ) 

LET NORM. CONSTANT = REMAI NOER/ I NTEGRAL 
FOR 1=1 TO PAIRS, DO 

LET RATE. FUNCT ION (I)= CUS TOME R (I ) *N0 RM. CONSTANT 
LOOP 

LET ‘I = INTEGRAL/CLOSING. TIME 

RESERVE S(*) AND U(«) AS VOLUME <-I*10 

RESERVE PLOT!'^) AS PAIRS 

LET RATE. NO = 1 

RETURN 

END 



ROUTINE ARRIVAL 
ADO 1 TO NO 
IF NO = 1 

LET S( 1 )=-LOG.E.F(PERCENT ) 

LET um =S(1 >/RATE.FUNCTI0N(1) 

IF U( 1) >FUNCTION( 1) 

LET RECORD = F UNO T I ON ( i » A T E . FUNC T I ON ( 1 ) 
GO TO TRANSFORM 
OTHERWISE 

ADD 1 TO PLOT(RATE.NO) 

RETURN W ITH U (1 ) 

ELSE 

LET S(NO)= -LOG.E.F(PERCENT) 

LET U(NO) = S( NO)/RATE.FUNCTION(RATE.NO) +U(N0-1) 

LET S(NO)=S(NO)i-S(NG-l ) 

HERE 

IF U(NO)>FUNCTION(RATE.NO) 

LET RECCRD= (FUNCT ION! RAT E. NO ) -FUNCT I ON ( RAT E. NO- 1 ) )* 
RATE.FUNCTION(RATE.NG) *■ RECORD 
• TRANSFORM' 

ADD 1 TO RATE. NO 
IF RATE. NO > PAIRS 

LET CASES = NO-lfNUM! 1) +NJM ( 2) <-NUM( 3) 
PRINT I LINE WITH CASES AND DAY THUS 
♦ ARRIVALS FOR DAY « 

RETURN W ITH U ( NO ) 

ELSE 

LET U( NOJ = ( S ( NO) -RECORD) /RATE. FUNCT ION! RATE.NO) + 
FUNCTION! RATE. NO-1 ) 

JUMP BACK 

ELSE 

IF UINO) >=CLOS ING.TIME 

LET CASES = NO-H-NUM! 1 ) i-NUM! 2 ) +NUM! 3 ) 

PRINT 1 LINE WITH CASES THUS 
ARRIVAL NO 
ALWAYS 

ADD 1 TO PLOTIRATE.NO) 

RETURN WITH U!NC) 

END 
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ROUTINE ACALC 

DEFINE I AS AN INTEGE R, SAVED VARIABLE 
IF TIME.V <= FUNCTION! 1) LET I = 1 REGARDLESS 

HERE 

IF TIME.OF.ARRIVAL(PATIENT) < FUNCT I ON ( I )-5 . 0 

RETURN 

ELSE 

IF TIME. OF. ARRIVAL! PATIENT) < FUNCT I ON ! I ) ♦■S . 0 
ADD 1 TO COUNTl ! I ) 

LET WAITING. TIME = TIME.V -T I M E . OF . A RR I V AL ! PAT I ENT ) 
ADD WAITING. TIME TO A. ME AN. OE L AY ! I ) 

ADD WAITING. TIME=:‘WAITING. TIME TO STATIlitI) 

IF WAITING. TIME > DEL AY . C R I T E R ION 

ADD 1.0 TO A. DELAY ! I ) 

ALWAYS 

RETURN 

ELSE 

ADD 1 TO I 
JUMP BACK 
END 



ROUTINE PCALC 

DEFINE I AS AN I NTE GER t SA VED VARIABLE 
IF TIME.V <= FUNCTION! 1) LET I = 1 REGARDLESS 

HERE 

IF TIME.OF.ARRIVAL!P>ATIENT) < FUNCTI ON ! I ) -5 . 0 

RETURN 

ELSE 

IF TIME.OF.ARRIVAL(PATIENT) < FUNCT I ON ( I )+5 . 0 
ADD 1 TO C0UNT2! I ) 

LET WAITING. TIME = TIME.V -T I ME .OF . ARR I V AL ! PAT I ENT ) 
ADO WAITING. TIME TO P . ME AN . D EL AY ! I) 

ADD WAITING. TIME-WAITING. TIME TO STAT!12,I) 

IF WAITING. TIME > DE LA Y . C R I TER I ON 

ADD 1. 0 TO P. DELAY! I ) 



ALWAYS 

RETURN 



ELSE 

ADD 1 TO I 
JUMP BACK 
END 
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PREAMBLE 

LAST COLUMN IS 60 



TEMPORARY ENTITIES 

EVERY AMOSIST MAY BELONG TO AN A VA I L . NEXT . AMOS AND HAS 
AN OFFICE 

EVERY PHYSICIAN MAY BELONG TO A NEXT .AVAIL . DOCTOR AND HAS 
A CUBICAL 

EVERY PATIENT MAY BELONG TO AN AMOS.CUEUE,A PHYS.QUEUEi 
A KEF.QUcUE, HAS A TYPE, A T I ME . OF . ARR I V A L AND A 
PH. ASSIST 

EVENT NOTICES INCLUDE END . OF . S I MUL AT I ON , CLOSE . CL I NIC , ■ 
OPEN. CLINIC, TRIAGE, UTIL IZAT ION AND CONSULT 
EVERY AMOS. APPOINT HAS A NUMBER 
EVERY PHYS. APPOINT HAS A LISTING 
EVERY CHG. IN. MEDICS HAS A SUBSCRIPT 
EVERY CHG. NUM. PHYS HAS AN INDEX 

THE SYSTEM OWNS A PHYS. QUEUE, AN AMOS. QUEUE, A 

NEXT. AVAIL. DOCTOR, AN AVA I L . NEXT. AMOS AND A 
REF. QUEUE 

DEFINE AVAI L. NEXT. AMOS , NEXT . AV AI L . DOCTOR AND AMOS. QUEUE 
AS FIFO SETS WITHOUT FF,FB,FA AND RL ROUTINES 
DEFINE PHYS. QUEUE AS A FIFO SET WITHOUT FF,FA AND RL 
ROUTINES 

DEFINE REF. QUEUE AS A SET RANKED BY LOW TI ME .OF . ARRI VAL 
WITHOUT RF AND RL ROUTINES 



PRIORITY 

-DEFINE 



DEFINE 
DEFI NE 

DEFINE 



DEFINE 

DEFINE 

DEFINE 

TALLY AV 
TALLY 
DEFI 
TALLY 

DEFI 
TALL Y 
TALLY 
DEFI 
DEFI 



ORDER IS END.OF .SI MULATION,CHG.IN.MEDICS , 
CHG.NUM.PHYS, AMOS. APPOINT , PHYS. APPO I NT, CONSULT, 
UTILIZATION, CLOSE. CLINIC AND TRIAGE 
NUMBER, LI STING, OFF ICE, CUBICAL , PA IRS, I NDEX , REF ERAL , 
SUBSCR I PT,NONREF , CONSOL, NONCONSCL, DAY ,NO,L AST. DAY, 
TYPE, VOLUME , R AT E . NO , NO. DUTY . CHGS , PH . A SS I ST AND 
NO. ROSTER. CHGS AS INTEGER VARIABLES 

NUM, NO. DOCTORS, NO. AMOSISTS AS I NTcGE R , 1 -DI M ARRAYS 
ALPHA, BET A, CLOSING. TIME , T IME . OF . ARR I V AL , RECORD , 
INTEGRAL AND DE LAY . CR I TE R I ON AS REAL VARIABLES 
S,U, TIME. TC, DOC TOR. T I ME. TC, FUNCTION, CUSTOMER, T I HE , 
AVE. ARRIVAL, PLOT, ENCOUNTER ANC RATE . FUNC T I UN AS 
1-DIM ARRAYS 

PERCENT AS A REAL FUNCTION 
ARRIVAL AS A REAL FUNCTION 
DPLTP AS A FORTRAN ROUTINE 



DEFI 
DEF I 



E. SERVICE. A AS THE MEAN OF A. SERVICE 
AV. SERVICE. P AS THE MEAN OF P. SERVICE 
NE A. SERVICE AND P. SERVICE AS REAL VARIABLES 
ACASE AS THE MEAN, SCASE AS THE STD. DEV AND MCASE 
THE MAXIMUM OF CASES 
NE CASES AS AN INTEGER VARIABLE 
MAXI AS TtiE MAXIMUM OF N. PHYS. QUEUE 
MEANl AS THE MEAN AND STD AS THE STD. DEV OF MAXD 
NE MAXD AS A REAL VARIABLE 

NE A. NO. DELAY, P. NO. DEL AY , A. DEL AY, P. DEL AY, A. UTIL, 
A.MEAN.DELAY,P .MEAN. DEL AY, P .UT I L, A. SYS. S I ZE, 

P. SYS. SIZE, A. AV. QUEUE AND P.AV. QUEUE AS 1-DIM 
ARRAYS 

NE MEAN AND STAT AS 2-DIM ARRAYS 
NE C0UNT1,C0UNT2 AS I NTEG ER , 1 -D I M ARRAYS 



AS 



BEFORE REMOVING FROM AMOS. QUEUE, CALL ACALC 
BEFORE REMOVING FROM PHYS. QUEUE, CALL PCALC 
END 



CASE 2 SIMULATION MODEL 
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MAIN 

CALL INPUT 
RELEASE INPUT 
AOO 1 TO DAY 
SKIP 1 OUTPUT LINE 
PRINT 1 LINE THUS 

SIMULATION BEGINS 
SCHEDULE AN OPEN. CLINIC NOW 
START NEW PAGE 
START SIMULATION 
END 



EVENT OPEN, 
DEFINE I 
DEFINE X 
RESERVE X 
SCHEDULE 
IF AVE.ARRI 
LET NLM( 1 
IF NUM(l) < 
LET NUM(3 
LET NUMd 
HERE 

IF AVE.ARRI 
LET NUMI2 
HERE 

FOR 1=1 TO 



FOR 1=1 TO 
CREATE 
IF T<= 



CLINIC SAVING THE EVENT NOTICE 
AS AN INTEGER VARIABLE 
AS A 1-DIM ARRAY 
( * ) AS 3 

A CLOSE. CLINIC AT CLOSING. TIME 
VAL(1)=0 JUMP AHEAD ELSE 
) =POISSON.F( AVE. ARRIVAL! 1) tl) 

= 0 JUMP AHEAD ELSE 
)=BINOMIAL. F( NUM( 1 ) , BETA, 1 ) 

)= (-NUM(3) ) +NUM( i ) 

VAL(2)=0 JUMP AHEAD ELSE 
J = P0ISS0N.F(AVE.ARRIVAL(2) ,1) 

3 DO 

LET X( I )=REAL.F(NUM( I ) ) 

LOOP 

NUM( 1 H-NUM( 2J +NUMI3) ,00 

LET L0AD=1.0/(X( ld-X(2)<-X(3) ) 



A PATIENT 
X( 1 )*LOAD 



LET T=PERCENT 

LET TYPE! PATIENT )=1 
LET X ( 1 ]=X( 1) - I 
FILE PATIENT IN AMOS. QUEUE 
CYCLE 



ELSE 

IF T<=(X(1} +X(2))*L0AD 



ELSE 



LET TYP£(PATIENT)=2 
LET X(2)=X(2)-1 
FILE PATIENT IN PHYS. QUEUE 
CYCLE 

LET TYPE(PATI£NT)=3 
LET X( 3) = X(3)-l 
FILE PATIENT IN AMOS. QUEUE 
LOOP 



SCHEDULE A TRIAGE AT ARRIVAL 

SCHEDULE A UTILUATION AT FUNCTION(l) 

SCHEDULE A CH G. I N. MED I CS NOW 

SCHEDULE A CHG.NUM.PHYS NOW 

LET INOEXICHG.NUM.PHYS ) = 1 

LET SUBSCRIPT (CHG. IN. MEOIuS ) = 1 

DESTROY OPEN. CLINIC 

RETURN 

END 
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EVENT CHG. IN. MEDICS SAVING THE EVENT NCTICE 
NORMALLY MODE IS INTEGER 
LET I = SUBSCRIPTICHG. IN. MEDICS) 

IF TIME.V=0.00 

FOR J = 1 TO NO.AMOSI STS{ 1) ,00 

CREATE AN AMOSIST 

LET CUBICAL (AMOSIST)=J 

FILE AMOSIST IN AVA I L . NEXT . AMOS 

LOOP 

GO TO APPOINTMENT 

ELSE 

LET M£DICS=NO.AMOSISTS( I ) -NO. AMOSI STS( I-l) 

IF MEDICS< 0 

FOR EACH AMOSIST OF A VA I L . NE XT . AMOS , DO 

IF CUBICAL (AMOSIST) >N0 . AMOS I ST S ( I ) 

REMOVE AMOSIST FROM AVAI L. NEXT . AMOS 
DESTROY AMOSIST 

ALWAYS 

LOOP 

JUMP AHEAD 
OTHERWISE 

FOR J = NO.AMOSISTS( I-D+l TO NO.AMOSI STS ( I ) ,00 

CREATE AN AMOSIST 

LET CUBICAL(AMOSI ST ) = J 

FILE AMOSIST IN AVAI L . NEXT . AMOS 

LOOP 

' APPCI NTMENT* 

IF THE AMOS. QUEUE IS NOT EMPTY 

LET AMOSIST = F . A VA I L . NEXT . AMOS 
REMOVE AMOSIST FROM AVA IL . N EXT . AMOS 
SCHEDULE AN AMOS. APPOINT NOW 
LET NUMBER (AMOS. APPCI NT) =CUBI CAL { AMOSIST) 
DESTROY AMOSIST 

ALWAYS 

HERE 

ADD 1 TO I 

LET SUriSCRIPT(CHG. IN. MEDICS) = I 
IF I <=NO.DUTY.CHGS 

THEN IF TIME. TCd XCLOSING. TIME 

SCHEDULE THIS CHG . I N. MED I C S AT TIME.TC(I) 

ALWAYS 

RETURN 

END 
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EVENT CHG.NUM.PHYS SAVING THE EVENT NOTICE 
NORMALLY MODE IS INTEGER 
LET I = INDEX(CHG.NOM.PHYS) 

IF TIME.V=0.00 

FOR J = i TO NO. DOCTORS! 1 ) »DO 

CREATE A PHYSICIAN 

LET OFFICE (PHYSICIAN) = J 

FILE PHYSICIAN IN NE X T . A VA I L . DOC TOR 

LOOP 

GO TO APPOINTMENT 

ELSE 

LET DOC = NO. DOCTORS! I ) -NO. DOCTORS ! I - 1 ) 

IF DOC <0 

FOR EACH PHYSICIAN OF NEX T . A V A I L . DOCT OR , DO 
IF OFF ICE! PHYSICIAN) > NO . DOC T OR S ( I ) 

REMOVE PHYSICIAN FROM NEXT . AVA IL . DOCTOR 
DESTROY PHYSICIAN 
ALWAYS 
LOOP 

JUMP AHEAD 
OTHERWISE 

FOR J = NO.DOCTORS! I-l) + L TO NO . DOCTOR S ! I ) t DO 

CREATE A PHYSICIAN 

LET OFFICE! PHYSICIAN) = J 

FILE PHYSICIAN IN NEX T . AVA I L . DOCTOR 

LOOP 

•APPOI NTMENT* 

IF THE PHYS. QUEUE IS NOT EMPTY 

LET PHYSICIAN = F .N EX T . AV A I L . DOC TOR 
REMOVE PHYSICIAN FROM NEXT . AV A1 L . DOCTOR 
SCHEDULE A PHYS. APPOINT NOW 

LET LISTING!PHYS.APPOINT)=OFF I CE ( PHY S I C I AN ) 
DESTROY PHYSICIAN 

ALWAYS 

HERE 

ADD 1 TO I 

LET INDEX!CHG.NUM.PHYS)= I 
IF I <= NO. ROSTER. CHGS 

THEN IF DOC TOR. TIME. TC! I ) < CLOSING. TIME 

SCHEDULE THIS CHG.NUM.PHYS AT DOC TOR . T IME . TC ! I ) 

ALWAYS 

RETURN 

END 
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EVENT TRIAGE SAVING THE EVENT NOTICE 
CREATE A PATIENT 

LET TIME. OF .ARRIVAL! PATIENT) = TIME.V 
LET T=PERCENT 
IF T <=ALPHA 

LET TYPc(PATIENT) =2 
FILE PATIENT IN PHYS. QUEUE 
JUMP AHEAD 
OTHERWISE 

IF T <= ALPHA 4- ( I .0-ALPHA J 1 ,0-BETAi 
LET TYPE! PATIENT)=1 
FILE PATIENT IN AMOS. QUEUE 
JUMP AHEAD 
ELSE 

LET TYPE!PATIENT)=3 
FILE PATIENT IN AMOS. QUEUE 
HERE 

IF TYPEIPATIENT) = 2 

THEN IF NEXT. AVAIL. DOCTOR IS NOT EMPTY 
LET PHYSICIAN = F . NEXT . A VA I L . DOCTOR 
REMOVE PHYSICIAN FROM NEXT .AVAI L . DOCTOR 
SCHEDULE A PHYS. APPOINT NOW 

LET L 1ST ING ! PHYS. APPOINT ) =OFF ICE! PHYSICIAN) 
DESTROY PHYSICIAN 
ALWAYS 

IF TYPEIPATIENT) = 1 OR TYPEIPATIENT) = 3 
THEN IF THE AVAI L. NEXT. AMOS IS NOT EMPTY 
LET AMOSIST = F. AV A I L . NEX T . AMOS 
REMOVE AMOSIST FROM AV A I L . NEXT . AMOS 
SCHEDULE AN AMOS. APPOINT NOW 
LET NUMBER ! AMOS. APPOINT) =CUB I CAL! AMOS I ST ) 
DESTROY AMOSIST 
ALWAYS 

LET NEXT = ARRIVAL 
IF NEXT < CLOSING. TIME 

SCHEDULE THIS TRIAGE AT NEXT 
RETURN 
OTHERWISE 
DESTROY TRIAGE 
RETURN 
END 
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EVENT AMOS. APPOINT SAVING THE EVENT NOTICE 
DEFINE I AS AN INTEGER VARIABLE 
LET I = NUMtsERIAMOS. APPOINT) 

IF I > NO. AM(jSISTS(SU3SCRIPT(CHG. IN. MEDICS)-!) 

DESTROY AMOS. APPOINT 
JUMP AHEAD 
ELSE 

IF THE AMOS. QUEUE IS EMPTY 

CREATE AN AMOSIST 

LET CUBICAL! AMOSI ST) =NUMBER (AMOS. APPOINT) 
FILE AMOSIST IN AV A I L . NEXT . AMOS 
DESTROY AMOS. APPOINT 
RETURN 

ELSE 

LET PATIENT = F. AMOS. QUEUE 

LET A.SERVICe= - ENCOUNTE R ( 1) *LOG. E . F ( PERCENT ) 

LET NEXT = A. SERVICE+T IME .V 
REMOVE PATIENT FROM AMOS. QUEUE 
IF TYPE(PATIENT)=I 

ADD I TO NONREF 
DESTROY PATIENT 

ALWAYS 

IF TYPE (PATIENT)=3 

ADD 1 TO REFERAL 

SCHEDULE A CONSULT AT NEXT 

LET T IME.OF .ARR I VAL( PATIENT ) =NEXT 

LET PH. ASSI ST( PATIENT)=NUMbER ( AMCS . APPO I NT ) 

DESTROY AMOS. APPOINT 

FILE PATIENT IN REF. QUEUE 

JUMP AHEAD 

OTHERWISE 

IF NEXT < CLOSING. TIME 

SCHEDULE THIS AMOS. APPOINT AT NEXT 

ALWAYS 

HERE 

IF THE AMOS. QUEUE IS NOT EMPTY 

FOR 1=1 TO N. AMOS. QUEUE 
WHILE AVAIL. NEXT. AMOS IS NOT EMPTY, CO 
LET AMOSIST = F . A V A I L . NEX T . AMOS 
REMOVE AMOSIST FROM AVA I L . N EX T . AMOS 
SCHEDULE AN AMOS. APPOINT NOW 
LET NUMBER (AMOS. A PPOI NT ) =CU Bi C AL ( AMOS I S T) 
DESTROY AMOSIST 
LOOP 

ALWAYS 
RETURN END 
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EVENT FHYS. APPOINT SAVING THE EVENT NOTICE 
DEFINE I AS AN INTEGER VARIABLE 
LET I = LISTING! PHYS. APPOINT) 

IF I > NO. DOCTORS! INDEX(CHG.NUM.PHYS)-l) 

DESTROY PHYS. APPOINT 
JUMP AHEAD 

ELSE 

IF THE PHYS. QUEUE IS EMPTY 

CREATE A PHYSICIAN 

LET OFFICE! PHYS I C I AN ) =L I S T I N G ! P HYS . APPO I NT) 
FILE PHYSICIAN IN NEXT .AVAIL . DOCTOR 
DESTROY PHYS. APPOINT 
RETURN 

ELSE 

LET PATIENT = F. PHYS. QUEUE 
IF TYPE!PAT1£NT)=E 

LET P. SERVICE = - ENCOUNT ER ( 3 ) *LOG . E . F ( PERCENT ) 

ADD 1 TO NONCONSOL 

ALWAYS 

IF TYPE(PATI£NTJ= 3 

LET P. SERVICE=- ENCOUNTER! 4) *L OG. E . F ! PERCENT ) 

ADO 1 TO CONSOL 

SCHEDULE AN AMOS. APPOINT AT P . SER V I C E *-T I ME . V 
LET NUMBER! AMOS. APPO I NT) = PH. ASS 1ST! PATIENT) 

ALWAYS 

REMOVE PATIENT FROM PHYS. QUEUE 
DESTROY PATIENT 

IF P. SERVICE + TIME.V < CLOSING. TIME 

SCHEDULE THIS PHYS. APPOINT AT P . SE RVI CE<-T I ME .V 

ALWAYS 

HERE 

IF THE PHYS. QUEUE IS NOT EMPTY 

FOR 1 = 1 TO N. PHYS. QUEUE 

WHILE NEXT. AVAIL. DOCTOR IS NOT EMPTY, DO 
LET PHYSICIAN = F . N EXT . A V A IL . OOC TO R 
REMOVE PHYSICIAN FROM NE XT .AVAI L . DOCTOR 
SCHEDULE A PHYS. APPOINT NOW 

LET LISTING! PHYS. APPOINT )=OFF ICE {PHYSICIAN) 
DESTROY PHYSICIAN 
LOOP 

ALWAYS 
RETURN END 
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EVENT UTILIZATION SAVING THE EVENT NOTICE 
DEFINE I AS INTEGERf SAVED VARIABLE 
IF TIHE.V=FUNCT ION( i ) LET I =1 REGARDLESS 
LET Y = NO. AMOSISTSI SU6SCRIPT(CHG. IN. MEDICS) -1 ) 

LET DUMMY = N. AMOS. QUEUE + Y - N . AVA I L . NEXT . A M OS 
ADD DUMMY TO A . S YS . S I Z E ( I ) 

ADO DUMMY * DUMMY TO STATIUI) 

ADD N. AMOS. QUEUE TO A . A V . QUE UE ( I ) 

ADD N. AMOS. QUEUE * N. AMOS. QUEUE TO STAT(3,I) 

LET DUMMY = 1.0- N . A V A I L .NEX T . AMOS/ Y 
ADD DUMMY TO A.UTIL(I) 

ADO DUMMY * DUMMY TO STAT(5»I) 

IF AVAI L. NEXT. AMOS IS NOT EMPTY 
ADD 1.0 TO A. NO. DELAY! I ) 

REGARDLESS 

LET Y = NO.DOCTORSlINDEX(CHG.NUM.PHYS)-l) 

LET DUMMY = N.PHYS. QUEUE Y - N. NEXT . AVAI L. DOCTOR 
ADD DUMMY TO P . SY S . S I Z E ( I ) 

ADD DUMMY « DUMMY TO STAT(2,I) 

ADD N.PHYS. QUEUE TO P . AV .QUEUE! I ) 

ADD N.PHYS. QUEUE - N.PHYS. QUEUE TO STAT!4,I) 

LET DUMMY = 1.0- N . NEX T . AV A1 L . DOCTOR/ Y 
ADD DUMMY TO P.UTIL! I ) 

ADD DUMMY « DUMMY TO STAT!6,I) 

IF NEXT. AVAIL. DOCTOR IS NOT EMPTY 
ADD 1.0 TO P. NO. DELAY! I ) 

REGARDLESS 

LET Y = FUNCTION! I ) 

IF Y <= CLOSING. TIME 
ADO 1 TO I 

SCHEDULE THIS UTILIZATION AT FUNCTION!!) 

RETURN 

ELSE 

DESTROY UTILIZATION 

RETURN 

END 
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EVENT CLOSE. CLINIC SAVING THE EVENT NOTICE 
NOf^MALLY MODE IS INTEGER 
IF THE REF. QUEUE IS NOT EMPTY 
PRINT 1 LINE WITH N. REF. QUEUE THUS 

AT CLOSING. TIME « PATIENTS ARE WAITING FOR CONSOLTATION 
RETURN ELSE 

PRINT 1 LINE WITH NUMi 1 ) , NUM ( 2 ) , NUM( 3 J THUS 
PRELOAOING FOR QUEUES I TO 3 IS « * * 

PRINT 1 DOUBLE LINE WITH REF ERAL , NONREF , CONSOL , NONCONS OL 
AND TIME.V THUS 

REFERALS * NONREF * CONSOLS * NONCON * CLOSING 
T X E 

FOR EACH AMOSIST OF AV A I L. NEXT . AM GS ♦ DO 

REMOVE AMOSIST FROM AVAI L . NEXT. AMOS 

DESTROY AMOSIST 

LOOP 

FOR EACH PHYSICIAN OF NEXT . AVAI L . DOC T OR , DO 

REMOVE PHYSICIAN FROM NE X T . A VA I L . DOCTOR 

DESTROY PHYSICIAN 

LOOP 

IF THE AMOS. QUEUE IS NOT EMPTY 

FOR EACH PATIENT IN AMOS. QUEUE, DO 

REMOVE PATIENT FROM AMOS. QUEUE 

DESTROY PATIENT 

ADD 1 TO A. PATIENTS. NOT. SEEN 

LOOP 

ALWAYS 

IF THE PHYS. QUEUE IS NOT EMPTY 

FCR EACH PATIENT IN PHYS. QUEUE, DO 

REMOVE PATIENT FROM PHYS. QUEUE 

DESTROY PATIENT 

ADD 1 TO P. PAT lENTS. NOT. SEEN 

LOOP 

ALWAYS 

PRINT 1 LINE WITH A . PA T I ENT S .NO T . S EEN AND 
P. PATIENTS. NOT. SEEN THUS 

LEFT IN AMOSIST CLINIC LEFT IN DOCTORS OFFICE 

PRINT 1 LINE WITH MAXI THUS 

THE MAX NUMdER OF PATIENTS IN THE PHYSICIAN QUEUE IS * 

SKIP 1 OUTPUT LINE 
LET MAXD = MAXI 
IF CAY =LAST.DAY 

PRINT 2 LINES WITH AV. SERVICE. P AND AVE . SERV ICE . A THUS 
AVERAGE PHYSICIAN SERVICE AVERAGE AMOSIST SERVICE 

PRINT 1 LINE WITH MEANl AND STD THUS 
THE AVE MAX FOR PHYS. QUEUE IS STD DEV 

SKIP 1 OUTPUT LINE 

PRINT 1 DOUBLE LINE WITH ACASE,SCASE AND MCASE THUS 
MEAN NO OF PATIENTS STD DEV MAXIMUM 

« PER CAY 

SCHEDULE AN END . OF .S I MUL AT I ON NOW 
RETURN 
ELSE 

LET TIME.V =0.00 LET DAY=DAY+1 LET N0=0 LET RAT£.NO=l 
LET REFEPAL=0 LET C0NS0L=0 LET NONCONSOL =0 LET NUNREF=0 
LET REC0RD=0 RESET TOTALS OF N. PHYS. QUEUE 
DESTROY CHG. IN. MEDICS DESTROY CHG. NUM. PHYS 

DESTROY CLOSE. CLINIC SCHEDULE AN OPEN. CLINIC NOW 

RETURN 

END 



84 



I 



I 



I 



( 



j 



EVENT END. OF. SIMULATION 

DEFINE MODCUR AS AN INTEGER VARIABLE 
PRINT 1 LINE WITH DAY THUS 
END. OF . SIMULATION DAY IS 

FOR J = 1 TO 0 ALSO FOR I = 1 TO PAIRS, DO 

IF J = 7 OR J = 6 LET STATU, II = MEAN! J, I) REGARDLESS 
LET MEANUtlJ = MEAN (J , I ) /L AS T . 0 AY 
LET STATUtl) = ST A T (J , I ) / L AS T . 0 AY-ME AN tJ , I) «ME AN (J , I) 
LET STATU, IJ = SORT . F ( STAT U , II ) 

LOOP 

FOR I = 1 TO PAIRS, DO 

IF COUNTim =0.0 JUMP AHEAD ELSE 
LET STAT(9,IJ = MEANI9, I ) 

LET MEAN(9,I) = HE AN( 9 , I ) /COUNT 1 ( I ) 

LET STAT ( 9, I )=SQRT.F( STAT( 9, I )/COUNTi( I ) - MEAN(9,D* 
MEAN (9, I J ) 

LET ME AN (II, I) = ME AN ( 1 1 , I ) /COUNT 1 ( I ) 

LET STATI 1 1 ,I J=SQRT.F(STAT( il,I J / COUNT 1 ( I ) -M EAN ( H , I ) 
MEANdl, I n 

HERE 

IF C0UNT2(IJ = 0.0 CYCLE ELSE 
LET STAT( 10,1 ) = MEAN! 10,1 ) 

LET M£AN(10,I) = MEANdO,!) /COUNTZdJ 
LET STAT( 10, I J=SQRT .F(STAT( 10, I ) /CCUNT2( I )-MEAN( 10, I } ={= 
MEAN! 10,1)) 

LET HEAN(12,I) = ME AN ( 1 2 , I ) /C OUNT 2 ( I ) 

LET STAT (12 , I ) =SQRT . F ( S TAT( 1 2 , I )/C0UNT2( I ) -M EAN (1 2 , I ) ^ 
ME AN ( 12,1 ) ) 

LOOP 

START NEW PAGE 
PRINT 1 LINE THUS 

AMOSIST SYSTEM 

SKIP 1 OUTPUT LINE 

PRINT 1 DOUBLE LINE WITH DEL AY . C RI TE R I ON THUS 

SYS SUE STD DEV AVE QUEUE STD DEV UTILIZ STD DEV 
P(0=0) STD DEV P(D> =!=.) AVE DEl.AY STD DEV NO TIME 
SKIP 1 OUTPUT LINE 
FOR I = 1 TO PAIRS, DO 

PRINT 1 DOUBLE LINE WITH ME AN (1 , I ) , ST AT ( 1 , 1 ) , ME AN ( 3 , I ) , 
ST AT (3, I) ,MEAN(5, I ),STAT(5, I ) ,MEAN(7, I ), ST AT (7, I) , 
KEAN(9, I ) ,MEAN( 1 1 , I ) , ST AT ( 1 1 , I ) , COUNT U I ) AND 
TIMEd) THUS 

:{c 3!c^3!s«3it :{s^3}:«3;< sjc ^ :«c 3j: 3lc 3jc ^ ' 

IF FRAC .F( 1/5 .0) =0.0 SKIP 1 OUTPUT LINE ALWAYS 
LOOP 

START NEW PAGE 
PRINT 1 LINE THUS 

PHYSICIAN SYSTEM 

SKIP 1 OUTPUT LINE 

PRINT 1 DOUBLE LINE WITH DEL AY. C RI TER I ON THUS 

SYS SIZE STD DEV AVE QUEUE STD DEV UTILIZ STD DEV 
P(C=0) STD DEV P(D> =i=.) AVE DELAY STD DEV NO TIME 
SKIP 1 OUTPUT LINE 
FOR I =1 TO PAIRS, DO 

PRINT 1 DOUBLE LINE WITH ME AN( 2 , I ) , STAT ( 2 , I ) , ME AN ( A , I ) , 
STAT( A, I ) , MEAM( o, I ) , STAT (6, I ) ,MEAN( 8 , I ) , STAT( 8,1), 
MEAN! lU, I ) ,MEAN(12, I ) , STAT( 12 ,9) ,CGUNT2( I ) AND 





TIME( I) THUS 
^ ^ 






❖ , >;c ^ 




^ ^ ❖ 


3}c ^ 


yic ^ 


.ir ^ « 


1 F 


FRAC.fi 1/5.0) 


=0.0 SKIP 


1 OUTPUT 


LINE ALWAYS 



LOOP 

START NEW PAGE 
PRINT 1 LINE THUS 

ACTUAL NUMBER OF ARRIVALS PER PERIOD 
SKIP 1 OUTPUT LINE 

FOR 1= 1 TO PAIRS, WRITE PLOT(I) AS D (15,4) 

FOR 1= 1 TO PAIRS-i, DO 

IF I = 1 LET PLOT ( I )=PLOT( 1 )* 60 . 0 / (FUNCT I ONI 1 )«LAST . DAY) 
ALWAYS 

LET PLOT! dl ) = PLOT (dl )« 60 . 0 /( (FUNCTION! I i- 1 ) - 
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FUNCTION( I ) }«LAST.DAY) 

LOOP 

CALL DPLTP(TIME1«) , PLOT < ^ , P A I RS , MCDCU R 
STOP END 



ROUTINE PERCENT 
LET X=RANDCM.F<2) 
RETURN UITH X 
END 



ROUTINE INPUT 

DEFINE I AS AN INTEGER VARIABLE 
DEFINE HOUR AND CLOCK AS 1-DIM ARRAYS 
DEFINE OPTION AS AN ALPHA, 1-DlM ARRAY 

RESERVE ENCOUNTER(=!') , AVE.ARRIVAL(>^) , NUM(«) AND OPTION!'^) 
AS 4 

LET ENCCUNTER(l) = 17.96 LET ENCOUNTER(2) = 23.83 
LET ENC0UNTER(3) = 13.55 LET ENC0UNTER(4) = 4.34 
LET VOLUME = 15o 

LET OPENING. TIME = 800. LET CLOSING. TIME = 2400. 

LET ALPHA = .55 LET BETA = .31 LET LAST. DAY = 1 

LET AVE. ARRIVAL! IJ = 7. LET AVE. ARRI VAL( 21 = 4. 

LET DELAY. CRI TERI0N=15.0 
HERE 

IF MODE IS NOT ALPHA 

PRINT 1 LINE WITH READ.V THUS 

EXPECTED OPTION WORD WHILE READING CARD NUMBER * 
STOP 
OTHERWISE 

IF EFIELD.F-SFIELD.F > 12 

PRINT 1 LINE WITH READ.V THUS 
WHILE READING CARO NUMBER * OPTION WORD WAS TOO LONG 
STOP 
OTHERWISE 

READ OPTION! 1 ) , OPT I ON! 2 ) , OPT I ON ! 3 ) 

IF OPTION! U = “ENCO" 

FOR I = 1 TO 4, DO 

READ ENCOUNTER! I ) 

LOOP 

JUMP BACK 

ELSE 

IF GPTI0N!1)= "VOLU" 

READ Y 

LET VOLUME= Y 
JUMP BACK 

ELSE 

IF OPTION! 1) = »OPEN" 

READ OPENING. TIME 
JUMP BACK 

ELSE 

IF CFTICNm =‘'CLOS» 

READ CLOSING. TIME 
JUMP BACK 

ELSE 

IF OPTION! 1) = ''PERC» 

READ ALPHA 
JUMP BACK 

ELSE 

IF CFTICN! 1) = ''REFE» 

READ BETA 
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JUMP BACK 

ELSE 

IF 0FTI0N(1}="ITER" 

READ Y 

LET LAST. DAY = Y 
JUMP BACK 

ELSE 

IF CPTI CN( 1 i =»PREL" 

READ AVE.ARRIVALC 1) , A VE . ARR I VAL ( 2 J 
JUMP BACK 

ELSE 

IF UPTION(I}="AMOS" 

READ Y 

LET NO.DUTY.CHGS = Y 

RESERVE TIME. TC tNO.AMOSI STS('-^) AMD CLOCK(«) AS 
NO.DUTY.CriGS 

FOR I = I TO NO.DUTY.CHGS ,DO 

READ CLOCK( I) fNO.AMOSI STS( li 

LOOP 

JUMP BACK 

ELSE 

IF GPTICNd ) = "PHYS" 

READ Y 

LET NO. ROSTER. CHGS = Y 

RESERVE DOCTOR. TIME. TC(** tNO.DCCTORS( ,HOUR(’i'} AS 
NO. ROSTER. CHGS 

FOR I = 1 TO NO. ROSTER. CHGS, DO 

READ HOUR{ I ) ,NO. DOCTORS! II 

LOOP 

JUMP BACK 

ELSE 

IF CFTICN{I)="ARRI" 

READ Y 

LET PAIRS = V 

RESERVE FUNCTION!*} ,R ATE . FUNCT I ON ! * i , 
CUSTOMER!*) tCOUNTI !*} ANDCOUNT2!*) AS PAIRS 
RESERVE TIME AS PAIRS 

RESERVE STAT!*,*} AND MEAN!*,*) AS 12 BY PAIRS 
FOR I = 1 TO PAIRS, DO 

READ TIME! I ), CUSTOMER! I ) 

LOOP 

JUMP BACK 

ELSE 

IF OPTION!!) ="DELa" 

READ DELAY. CRITERION 
JUMP BACK 

ELSE 

IF CPTICN!!) -.= "END." 

PRINT 1 DOUBLE LINE WITH OPTION!!) THUS 
**** IS NOT THE PkOPER BEGINNING TO AN OPTION WORD. SIMULA! 
ION ENDED 

STOP 

OTHERWISE 

PRINT 1 DOUBLE LINE WITH ENCOUNTER!!), ENCOUNT ER ! 2 ) , 
£NCOUNTER!3) AND £NC0UNTER!4) THUS 
ENCOUNTER!!) ***.** ENCOUNTER!2) ***.** ENCOUNTER!3) ***. 
** ENCOUNTER! 4) *.** 

PRINT 1 DOUBLE LINE WITH VOLUME , 0 PEN I NG . T I ME , CLOS ING .T I ME , 
ALPHA, BETA AND LAST. DAY THUS 

VOLUME * OPENING TIME * CLOSING TIME * PERCENT TO 

DOCTOR .*** PERCENT REFERRED .*** NO OF DAYS * 

SKIP I OUTPUT LINE 
FOR I = I TO 2 , DO 

PRINT ! LINE WITH I , AVE . ARR I VAL! I ) THUS 
AVE PRELOADING OF QUEUE * IS * 

LOOP 

SKIP I OUTPUT LINE 
IF NO.DUTY.CHGS <= 0 

LET NO.DUTY.CHGS = 5 

RESERVE T IME.TC!*) ,NO.AMOSISTS!*) AND CLOCK!*) AS 5 
LET NO.AMOSISTS! 1 ) =5 LET NO . AMOS I ST S ! 2) =3 LET 

NO.AMOSI STS!3) =5 LET NO . AMOS IS T S ! 4 ) =7 LET 



87 



NO. AMOSI STS( 5 J =2 LET CLOCK(l) =800. LET CLOCK ( 2 J = 1 1 00. 
LET CLOCK ( 3) =1300. LET CLOCK ( 41 = 1 530 . LET CLOCK ( 5 1 = 1 6 :sO. 
REGARDLESS 

IF NC.kCSTER.CHGS < = 0 

LET NO. ROSTER. CHGS=5 

RESERVE DOCTOR. TIME. TC(=f‘) , NO. COCTCRS(=!=), HOUR!*) AS 5 
LET NO. DOCTORS! 1 J =2 LET NO . DOCT 0 RS ( ^} = 1 LET 
NO. DOCTORS! 3) =2 LET NO . DOC TOR S ! 4 i =3 LET NO. DOCTORS ! 5 1 =1 
LET HOUR!11=800. LET HOUR ! 2 ) =1 1 00 . LET HOUR ! 3 ) =1 2 00 . 
LET HOUR! 4) =1530. LET HOUR ! 5 ) = 1 63 0. 

REGARDLESS 

FOR I = 1 TO NO.DUTY.CHGStDO 

LET Y = CLOCK! IJ-OPENING. TIME 
LET Y=Y/100.0 

LET TIME.TL!!) = TRUNC . F ! Y J * 60 . 0 + FRAC .F ! Y) 100 . 0 
PRINT 1 LINE WITH NO. AMOS 1ST S ! I) , CLOCK! I I THUS 
NUMBER OF AMOSISTS ON DUTY *>!'=;< TIME SHIFT STARTS * 

LOOP 

SKIP 1 OUTPUT LINE 

FOR I = 1 TO NO. ROSTER .CHGS, DO 

LET Y = HOUR! n -OPENING. TIME 
LET Y= Y/100.0 

LET DOCTOR. TIME. TCm=TRUNC.F!Y)=i'60.0 +FRAC. F ! Y ) =!= 
100.0 

PRINT 1 LINE WITH NO .DOC TORS ! II ,HOUR ! I ) THUS 
NUMBER CF DOCTORS ON DUTY TIME SHIFT STARTS 

LOOP 

SKIP 1 OUTPUT LINE 
IF PAIRS<= 0 

LET PAIRS=32 

RESERVE FUNCTION!*) t RATE . FUNC T ION !* ) , CUSTOMER!*) , 

COUNTl !*) ,C0UNT2! *) AND TIME!*) AS 32 
RESERVE STAT!*,*) AND MEAN!*,*) AS lA BY 32 
LET TIME!1)=810. LET TIME!2)=820. LET TIM£!3)=840. LET 
TIME!4)=850. LET TIME!5)=900. LET TIME!6)=910. LET 
TIME!7)=920. LET TIME!8)=930. LET TIME!9)=945. LET 
TIME! 10) =1000. LET T I ME ! 1 1 ) = 1 0 1 5 . LET T I ME ! 1 2 ) = 1030 . 

LET TIME ! 1j) = 1045. LET T IHE ! 14) =1 100. LET TIME!15) = 
1115. LET TIME! 16) =1130. LET TI ME ! 1 7 ) =1 1 45. LET 
TIME! 18) = 1200. LET TI ME ! 19 ) = 1 21 5 . LET T I ME ! 2 0 ) =1245 . 
LET TIME ! 21 ) = 1300. LET T IM E ! 2 2) = 1 530. LET TIME!23) = 
1400. LET TIME!24) = 1430. LET T I ME ! 25 ) =1 50 0 . LET 
TIME! 26) =1530. LET TI ME ! 27 ) = 1 630 . LET TI ME ! 28 ) = 1 730 . 

LET TIME! 29) = 1930. LET T I ME 1 30) =2 130. LET TIME!31) = 
2330. LET TIMEI32 )=2400. 

LET CUSTOMER! 1)=22. 5 LET CUSTQME R ! 2) = 25. 5 LET CUST0MER!5)= 
24.5 LET CUSTOMER! 4)=23. 5 LET CUSTOMER ! 5 )=22. 5 LET 
CUST0MER!6) = 21.5 LET C US TOMER ! 7 ) = 20. 5 LET OUST OME R ! 8 ) = 1 9. 6 
LET CUSTCM£R!9)=18.6 LET CU STOME R ! 10 ) = 1 7. 6 LET CUSTOMER! 11) 
=16.0 LET CUSTOMER! 12) =14. 7 LET C U ST OMER ! 1 3 ) = 1 5 . 3 5 LET 
CUSTCMER!14)=12.05 LET C USTOMER ! 1 5 ) = 1 0 . 95 LET CUSTOMER ! lo) = 
9.9 LET CUST0MER!17)=8.9 LE T C USTOMER ! 1 3 ) =7. 9 LET 
CUSTOMER! 19) =7. LET CU ST OMER ! 20 ) =6 . 2 LET CUSTOMER ! 2 1 ) = 7 . 2 
LET CUSTOMER! 22)=13. 6 LET CU STOME R! 25 ) = 12. 5 LET CUSTGMER!24) 
=11.1 LET CUSTOMER ! 25) =9.85 LET C USTOMER ! 26) =8 . 8 LET 
CUSTCMER!27)=7. 55 LET C USTOM ER ! 28 ) =6 . 3 LET C U S TOMER ! 29 ) = 5. 8 
LET CUST0MER!30) =4.4 LET C USTOME R ! 51 ) = 2 . o LET CUSTOMcR! 32 ) = 
.01 

REGARDLESS 

FOR 1= 1 TO PAIRS, DO 

LET Y =!TIME! I )-OPENING.TIME)/100.0 

LET FUNCTION! I )=TRUNC.F!Y)*60.0 + FRAC.F!Y) *100.0 
PRINT 1 LINE WITH I , C USTOME R ! I ) , T I ME ! I ) THUS 
PERIOD NUMBER * AVE ARRIVAL RATE *.^<‘=!= AT TIME * 

LOOP 

LET Y = ICLOSING. TIME-OPENING. TIMEl/lOO.O 

LET CLOSING. TIME=TRUNC .F !Y)*60.0 F R AC . F ! Y ) * 1 00 . 0 

LET A.SY5.SIZE! *) = MEAN!1,*) LET P . SYS . S I ZE ! * ) = MEAN!2,*) 

LET A. AV.QUEUE! * ) = MEAN!3,*) LET P . A V .CUEUE ! * ) = MEAN!4,*) 

LET A. util!*) = MEAN!5,*) LET P.UTIL!*) = M£AN!6,*) 

LET A.NO.UELAY !* ) = MEAN!7,*) LET P . NO . DE L AY ! * ) = MEAN!8,*) 
LET A. DELAY!*) = MEAN! 9,*) LET P. DELAY!*) = MEAN!10,*) 
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LET A,MEAN.DELAY(*} = MEANdl,*) LET P . ME AN . D EL A Y ( * J = 

MEAN ( 12 »* ) 

IF PAIRS=1 

LET INT£GRAL=CLOSING.TIME«CUSTOMER( I 1 
JUMP AHEAD 

ELSE 

LET INTEGRAL = FUNCT I ON ( 1 } ^CUSTOM £R ( 1 ) 

FOR I = 2 TO PAIRS WHILE FUNCTION! I) < CLOS ING .T I ME t DO 

LET INTEGRAL = i NT EGRAL *- ( FUNCT I CN( I J -F UNC T I ON ( I- i ) ) 
^CUSTOMER ( I ) 

LOOP 

LET INTEGRAL = I NT EGRAL+ ( CLOSI NG . T I ME-FUNCTI ON ( I - 1 )} 
^•CUSTOMER ( I ) 

HERE 

LET REMAINDER=VOLUME-AVE. ARRIVAL! 1 )-AVE. ARRIVAL! 2) 

LET NORM. CONSTANT = REMAINDER/ I NTEGRAL 
FOR 1=1 TO PAIRS, DO 

LET RATE. FUNCTION (I)= CUSTOMER! I } ’{'NORM. CONSTANT 
LOOP 

LET I = INTEGRAL/CLOSING. TIME 

RESERVE Sl’i'l AND U!’!') AS VOLUME ^dlO 

RESERVE PLQT(^) AS PAIRS 

LET RATE. NO = 1 

RETURN 

END 



ROUTINE ARRIVAL 
ADD 1 TO NO 
IF NO = 1 

LET S!l)=-LOG.E.F(PERCENT) 

LET Um =S m/RATE. FUNCTION! i) 

IF U! I) >FUNCTION(l ) 

LET RECORD = FUNCT I ON! 1 ) *RATE . FUNCT I ON (1 ) 
GO TO TRANSFORM 
OTHERWISE 

ADD 1 TO PLOT !RATE.NOJ 
RETURN WI TH U ! 1 ) 

ELSE 

LET S!NO)= -LOG. E.F! percent; 

LET U(NO) = S! NO}/RATE. FUNCTIONIRAT E.NO) +U1NO-1) 

LET S!NO)=S!NO)+S!NO-l) 

HERE 

IF U!NO}>FUNCTION!RATE.NO) 

LET RECORD=!FUNCTION(RATE.NO)-FUNCTI ON!RATE.NO-i ) )* 

rate.functionirate.no; + record 

•TRANSFORM* 

ADD 1 TO RATE.no 
IF RATE. NO > PAIRS 

LET CASES = NO-lt-NUM! 1) +NUM! 2; +NUM( 3 ; 
PRINT 1 LINE WITH CASES AND DAY THUS 
* ARRIVALS FOR DAY - 

RETURN WITH U!NO; 

ELSE 

LET U1N0;=(S! M0;-REC0RD; /RATE.FUNCTION(RATE.NO;+ 
FUNCTION! RATE .NO-I ; 

JUMP BACK 

ELSE 

IF U!NO; >=CLCSING.TIME 

LET CASES = NO-dNUM! 1 ; <-NUM! 2I+NUM! 3 ; 

PRINT 1 LINE WITH CASES THUS 
ARRIVAL NO 
ALWAYS 

ADO 1 TO PLOTIRATE.NO; 

RETURN WITH U!NO; 

END 
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ACALC 

E I AS AN INTEGER, SAVED VARIABLE 

.V <= FUNCTION! 1) LET I = I REGARDLESS 



ROUTINE 
DEF IN 
IF TIME 
HERE 

IF TIME.OF.ARRIVAL(PATIENT) < FU NC T I ON ( I ) -5 . 0 



ELSE 
IF TIME 



RETURN 



ELSE 
ADD I 
JUMP 
END 



.OF .ARRIVAL! PATIENT ) < FUNCT I ON ! I H-5 . 0 
ADO 1 TO COUNTl ! I ) 

LET WAITING. TIME = TIME.V -T I ME .OF . AR R I V AL ! PAT I ENT ) 
ADO WAITING. TIME TO A. ME AN. DE L AY ! I ) 

ADD WAITING. TIME*WAITING. TIME TO STAT!11,I) 

IF WAITING. TIME > DEL AY . CR I T ER I ON 

ADD I .0 TO A. DELAY! I) 

ALWAYS 

RETURN 

TO I 
BACK 



ROUTINE 
DEF IN 
IF TYPE 

ELSE 
IF TIME. 
HERE 
IF TIME 

ELSE 
IF TIME 



PCALC 

E I AS AN INTEGER, SAVED 
(PATIENT) = 3 
RETURN 



VARIABLE 



V <= FUNCT ION! 1 ) 



LET I = 1 REGARDLESS 



OF.ARRIVAL!PATIENT) < FUNCTI ON ( I ) -5.0 



RETURN 



ELSE 

ADO 

JUMP 

END 



0F.ARRIVAL!PATIENT) < FUNCTION!! )+5.0 
ADD 1 TO C0UNT2 !I ) 

LET WAITING. TIME = TIME.V -T I ME . OF . AR R I V AL ! PA T I EN T ) 
ADO WAITING. TIME TO P. MEAN. DELAY ! I ) 

ADO WAITING. T IME*WAITING. TIME TO STAT!12,I) 

IF WAITING. TIME > DE LAY . CR I T ER I ON 

ADO 1.0 TO P. DELAY! 1) 

ALWAYS 

RETURN 

TO I 
BACK 



EVENT CONSULT SAVING THE EVENT NOTICE 
NORMALLY MODE IS INTEGER 
DESTROY CONSULT 
IF TIME.V > CLOSING. TIME 

SCHEDULE A CLOSE. CLINIC NOW 

REGARDLESS 

LET PATIENT = F. REF. QUEUE 
REMOVE PATIENT FROM REF. QUEUE 
IF THE NEXT. AVAIL. DOCTOR IS NOT EMPTY 
FILE PATIENT IN PHYS.CUEUE 
SCHEDULE A PHYS. APPOINT NOW 
LET PHYSICIAN = F . N5X T . AV A I L . DOCT 0 R 
REMOVE PHYSICIAN FROM NEX T . A V A I L . UOCT OR 
LET LI STING! PHYS. APPOINT) = OFF I C E ! PH YS I C I AN ) 
RETURN 
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F,PHYS. QUEUE 



OTHERWI SE 
LET KIND = 

HERE 

IF KIND -^= 0 

IF TYPE(KIND) -. = 3 

FILE PATIENT BEFORE KIND IN PHYS. QUEUE 
RETURN 
OTHERWISE 

LET KIND = S. PHYS. QUEUE! KINDI 
JUMP BACK 

OTHERWISE 

FILE PATIENT IN PHYS. QUEUE 
RETURN END 
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9 






ENCOUNTER 

17.96 23.85 13.55 4.34 

VOLUME 

156 

OPENING. TIME 
0800 

CLOSING. TIME 
2400 

PERCENTAGE 

.55 

REFERRALS 

.31 

I TEKATIONS 
200 

PRELOAOING 

7.0 4.0 
AMOSIST.SCH 

5 0800 . 5 . 1100 . 3 . 1300 . 5 . 1530 . 7 . 1630 . 2 . 

PHYSICIANS 

5 0800 . 2 . 1100 . 1 . 1200 . 2 . 1530 . 3 . 1630 . 1 . 

ARRIVAL. FCN 

32 810 . 22.5 820 . 23.5 840 . 24.5 850 . 23.5 930 . 22.5 910 . 
21.5 920 . 20.5 930 . 19.6 945 . 18.6 1000 . 17.6 1015 . 16.6 
1030 . 14.7 1045 . 13.35 1100 . 12.05 1115 . 10.95 1130 . 9.9 
1145 . 8.9 1200 . 7.9 1215 . 7.0 1245 . 6.2 1300 . 7.2 1330 13.6 

1400 . 12.5 1430 . 11.1 1500 . 9.85 1530 . 8.8 1630 . 7.55 1730 . 

6.3 1930 . 5.8 2130 . 4.4 2330 . 2.8 2400 . .01 
DELAY. VALUE 

15.0 

END. OF. DATA 



SAMPLE INPUT 
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